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Abstract 



Path integration over Eudidean chiral fermions is replaced by the quantum mechanics of an 



Q_i auxiliary system of non-interacting fermions. Our construction avoids the no-go theorem and 

^ faithfully maintains all the known important features of chiral fermions, including the violation of 

some perturbative conservation laws by gauge field configurations of non-trivial topology. 
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1. Introduction 

A surge of interest [1-6] in the problem of constructing lattice chiral gauge theories has been 
generated recently by a proposal of David Kaplan [1]. The ultimate goal of this activity is to 
make a study of chiral gauge theories outside perturbation theory possible. While the simulation 
of pure gauge theories on the computer is conceptually easy, the inclusion of chiral fermions has 
been a theoretical stumbling block for years [7]. The nature of the stumbling block is essentially 
kinematical in that it only involves the fermion dynamics in a fixed gauge background, and this 
problem only requires to make sense of a quadratic (also infinite and Grassmann) path integral. The 
fermion dynamics results in an effective action (the logarithm of the chiral determinant) which is a 
complex functional of the gauge field background. A correct local lattice regularization of the chiral 
determinant should satisfy two main requirements: Firstly, it should reproduce perturbation theory 
(for small and smooth gauge fields). Secondly, it should reproduce instanton physics correctly (for 
large and smooth gauge fields that carry non-zero topological charge). In the present paper we 
present such a regularization. Having overcome the kinematical problem we hope to render the 
insufficiently understood dynamics of chiral matter coupled to nonabelian gauge fields amenable to 
non-perturbative study. There are further obstacles on the road to efficient numerical simulation; 
we hope to convince the reader that the seemingly insurmountable kinematical obstruction has been 
removed and therefore it is worth investing efforts into dealing with the more technical problems 
still ahead. 

Our regularization is based on a formal representation of the chiral determinant by the overlap 
of two states in an auxiliary Hilbert space. The two states are the many body ground states of 
two Hamiltonians describing non-interacting fermions. After regularization, in a finite Euclidean 
volume, the states become finite vectors with components labelled by lattice sites, spinorial and 
internal symmetry indices. The overall phases of the states are defined in a way that reflects several 
discrete symmetries of the continuum target theory. In [4], the overlap representation was arrived 
at in a rather convoluted manner: Kaplan[l], started from a Dirac fermion on a 2c? + 1 dimensional 
lattice with a mass term inhomogeneous in one direction creating a transverse domain wall. The 
Dirac operator has a zero mode with a definite chirality bound to the domain wall. Kaplan's 
proposal was to couple this zero mode to a gauge field in order to construct a 2d dimensional chiral 
gauge theory. At the same time there was an independent proposal by Frolov and Slavnov [8] to 
perturbatively regulate an anomaly free chiral gauge theory by using an infinite number of Pauli 
Villars regulators. It was shown in [3] that the above two proposals can be seen as related once 
it is recognized that the abstract central idea is to describe a single chiral fermion by an infinite 
flavor tower of Dirac fermions. The Dirac mass matrix is representable by a finite operator in 
infinite flavor space whose analytical index is equal to unity times a prefactor of the order of the 
ultra-violet cutoff. The infinity of the number of fermions can be rigorously controlled and the 
result is the overlap formula for the chiral determinant [4]. 

The overlap formula can be motivated directly, without recourse to an infinite number of 
fermions. In this paper we construct, in complete detail, chiral gauge theories on a lattice starting 
with the overlap formalism. Since there is no reference to an infinite number of Fermi fields, 
one does not have to worry about the apparent need to tackle the seemingly impossible task of 
incorporating them exactly. Neither does one stumble into any of the many pitfalls arising out 
of truncating them and spoiling the important property of the mass matrix, namely, the non-zero 
index. The direct motivation of the overlap formulation is elegant, natural and convincing. From 



this vantage point, the infinite number of fermions appears as an after thought, a possible new 
interpretation of the overlap. Still, we feel that it is necessary to keep both view points in mind 
for a complete understanding of our construction. For example, if we insist on a description of our 
regularization in a Hamiltonian (rather than Euclidean action) framework, we need locality in the 
physical time direction and are forced to think in terms of an infinite number of fermions. 

In this paper we focus on the lattice regularization of the fermionic part of the chiral gauge 
theory. The gauge part can be regulated on the lattice in standard manner. There are three parts 
to this paper. In the first part (sections 2-6), we deal with the overlap formula and the formulae for 
fermion expectation values. In section 2, we present the connection between the chiral determinant 
and the overlap formula at a formal level. The formal overlap is easier to regulate than the chiral 
determinant. The latter vanishes for the large class of gauge fields carrying non-zero topological 
charge. In section 3 we show how the overlap formula can easily reproduce this important property. 
The phase of the chiral determinant contains all the parity breaking effects induced by the chiral 
fermions and, in particular, possible anomalous breaking of gauge invariance. This role is taken up 
by the phase of the overlap which is sensitive to the phases of the two participating states. The real 
part of the overlap does not depend on these phases and its definition is naturally gauge invariant. 
Gauge invariance is broken by employing a Wigner-Brilouin phase choice for the many body ground 
states. In section 4, we show that this phase choice preserves several important formal properties 
of the imaginary part of the continuum induced action. Having defined the overlap, we proceed to 
construct representations of fermion expectation values in section 5. In particular, we also show 
how to define fermion expectation values in an anomaly free gauge theory when the background 
gauge field carries a non-zero topological charge. In section 6 we describe how to compute the 
chiral determinant and the fermion expectation values in perturbation theory. 

Having discussed the overlap formula and the fermion expectation values, we proceed to con- 
struct a lattice regularization of these quantities in the second part (section 7-9) of the paper. In 
section 7 we show that the two many body Hamiltonians whose many body ground states com- 
prise the overlap can be regulated in a straightforward manner. The regularization so obtained is 
simpler to understand and easier to work with than that obtained in [4], but both regularizations 
are conceptually equivalent. Section 8 deals with the vanishing of the overlap in non-trivial gauge 
backgrounds. An important byproduct is a clean definition of the topological charge for lattice 
gauge fields, this time based directly on the reaction of fermions. Having constructed the needed 
ingredients, regularized anomaly free theories are put together in section 9. Vector theories with 
arbitrary numbers of flavors are also similarly constructed. The fermionic propagators in the vec- 
tor theories satisfy the basic properties needed to prove the well known mass inequalities involving 
mesons and baryons. The proof is now in a rigorously regularized theory and holds for an arbitrary 
number of flavors, both for massless and for massive fermions. 

The third part (section 10-11) deals with some tests in two and four dimensions. In section 
10 we focus on two dimensional abelian models. We compute the overlap for perturbative gauge 
fields and show that both the magnitude and the phase agrees with the continuum result. The 
overlap reproduces the correct topological charge when the background gauge field carries one. We 
study the overlap in a gauge background made up of a localized instanton and a well separated 
localized anti-instanton. The quantitative behavior of the overlap as a function of the instanton 
anti-instanton separation is as expected from continuum considerations. The fermionic two-point 
function also has the expected continuum behavior. By explicitly looking at a background with a 
single localized instanton, we extract the zero mode and show that it has the right shape. Further, 



we study the axial anomaly in a vector theory in the presence of a background gauge field that carries 
non-zero topological charge. We find that the continuum anomaly still is reproduced, indicating 
that the definition of the overlap is also valid for non-perturbative gauge fields. All results are in 
complete agreement with the exact results for the Schwinger model and it is argued that there is a 
good chance this formalism will work if applied to a chiral rather than vector gauge theory in two 
dimensions. Four dimensional gauge theories are discussed in section 11. Since the computational 
effort is larger in four dimensions, fewer tests have been carried out so far. We show that the 
perturbative anomaly is reproduced and the overlap correctly accounts for the number of zero 
modes both in abelian and non-abelian backgrounds. The results of these important tests stongly 
indicate that this formalism should also work for chiral gauge theories in four dimensions. We end 
this paper with a discussion in section 12 of the expected behavior when the dynamics of the gauge 
fields are included. 



2. Chiral determinant as an overlap of Slater determinants 

The massless Euclidean Dirac operator, 

B = j^{d^ + iA^{x)), (2.1) 



in the chiral basis. 



a^\ fl 



is of the form 



where 



In 4-D 



and in 2-D 



C = a^{d^ + iA^{x)). (2.4) 

(Tl = 1, (72 = i. (2.6) 



A^(a;) is a finite hermitian matrix for each x and /x and this includes as special cases all compact 
gauge groups with the fermion being in any representation. 

In a vector theory, integration of the fermions in a fixed background gauge field yields det D 
which by virtue of (2.3) is 

det D = det C det C^ = det CC^ = det C^C = | det Cf (2.7) 

This shows that the vector determinant is a product of two chiral determinants of opposite chirality. 
D, CC^, C^C are operators that map vectors from one space into itself and therefore one can think 
of their determinants as products of eigenvalues and proceed to regularize these products. The same 



is not true of the determinants of C and C^^ individually. C is the chiral Dirac operator and is a 
mapping from Vl to Vr where elements of Vl are decreed to transform as {\,0) under the Euclidean 
Lorentz group and those of Vr are decreed to transform as (0,^) under the same Lorentz group. 
As such det C cannot be thought of as a product of eigenvalues. 

Let {ui} and {vj} be orthonormal bases of Vr and Vl respectively. These bases were picked 
without any reference to the gauge potential A. A reasonable definition, compatible with (2.7), of 
det C would be 

det C = det < Ui Cvj > . (2.8) 

ij 

This formal definition is also compatible with the known rules for Feynman graphs because if 
we expanded C around Cq = (7^5^ and chose the bases as plane waves with ordinary spinorial 
structure we indeed recover unrenormalized perturbation theory. Explicitly, up to normalization 
factors, (2.8) is an overlap of two many body wave functions*, i.e., 

det C = [[d^] [det(u,(^,))l * det((Ci;,(^,)). (2.9) 

The two many body wave functions in (2.9) are written as Slater determinants and the variables ^ 
denote all the continuous and discrete arguments of the single particle wave functions. 

We now proceed to describe one approach to obtain the two many body wavefunctions whose 
overlap is det C and where Lorentz invariance is maintained. The idea is to build a larger vector 
space, V = Vr ® Vl, and two hermitian operators in it that can be viewed as Hamiltonians, H : 

H±=(^7 ^^^=l5[±m + j,{d, + iA,{x))]; m>0 (2.10) 

H are Dirac Hamiltonians in [d+ 1) dimensions in the presence of a static (from the point of view 
of the (c? + 1) dimensional world) gauge field whose time component vanishes. Suppose now that 
the mass squared, m^, is much larger than ||CC^|| (at this point we are still at the formal level and 
the fact that ||CCl^|| = oo is ignored). The Dirac sea for H"*" (H~) is given by a Slater determinant 
made out of one particle wave functions with large bottom (top) components. This would be the 
ground state for a system of identical fermions (in the grand canonical ensemble and with zero 
chemical potential) all subjected to the single particle Hamiltonian H"*" (H~). If m is very large 
the many body ground states for H"*" and for H~ will be orthogonal. We cannot therefore ignore C 
entirely. Taking into account Fermi statistics we conclude that the overlap between the two ground 
states is linear in det C to leading order in - — 2~^. Since both ground states are Lorentz invariant 
objects we have a Lorentz invariant formula. Here we still refer to the four dimensional Lorentz 
group as the "Lorentz Group". 

We now derive the explicit connection between the overlap of the two many body states of H 



* This is equivalent to the mathematical definition [9] which regards det C G (det y2^)*®det Vr as 
the natural element corresponding to C G V^i^Vr. det ij(^Ui(^^j)) G det Vr and detij(^Vi(^^j)) G det Vl 
are two single element bases while the map det C between the two one dimensional spaces acts by 
(det C) detij(i;i(^j)) = det(Ci;i(^j)). (2.9) computes the (single) matrix element of det C relative 
to the two bases. 



and det C. The single particle eigenstates are given by 



yPj 



Y,ii^{xai,y(3j;A)i;f,^{y(3j;A) = \f,^i;f,^{xai; A); Ai± < 



(2.11) 



yPj 



a,P and i,j denote spin and color indices. Let |X± >a denote the many body states obtained by 

K 

(2.12) 



filling all the A^ states. We write 



i^k^ 



UK 
Vk 



and (2.10) implies 



Cvk = (A£.^ T rn)uK C^uk = {\ff ± m)vK 

Operating with C and C^ on the last two equations we get: 

C^CvK = ( [A^±] ' - m')vK CC V = ( [A^.^] ' - m')uK 



(2.13) 



(2.14) 



The eigenvalues \j^ of H^ are simply related to the eigenvalues of C^^C and CC^^ (the spectra 
of these two operators are identical except for possible zero modes). The manifold of all C's is 
parameterized by the set of all hermitian matrix valued vector potentials with no restrictions. As 
such, it is important that we write expressions for the eigenfunctions that are valid for all possible 
vector potentials, including the cases where C^^C and/or CC^^ have zero modes (in these cases 
H has some eigenvalues equal to m in absolute value). We note that \j^ < since the [d + 1) 
dimensional Dirac problems we are interested in always have a nonvanishing gap in the fermion 
spectrum. For H"*" the single particle wave functions of interest are 



V' 



K 



2^xx 



Cir^^x 



Vk 



2^xx 



2^xx 



^KK' =2 < ""K 



-Vk 

etc 



Cvk 



(VCtC + m2 +m)2 



Vk' > 



C 



■-i/ctC+^-m 

Vk 



Vk 



(2.15) 



For H the single particle wave functions of interest are 



i^k- 



2N 



UK 



-UK 



KK 



2N 



KK 



2N 



KK 



N 



KK' 



-\ < UK 



cct 



UK 
^ C^UK 



(VCCt + m2 +m)2 
5 



Ct- 



Uk 

1 



VcCt + r 



-UK 



UK' > 



(2.16) 



Note that the matrices N'^ and N are diagonal. 

The many body states |X± >a are Slater determinants of the single particle wave functions 
Tpj^ and the overlap 

A<L- \L+ >A = det (< V'^-IV'^t >) 

KK' 

1 , , 1 (2-17) 

= . = det < uk , ^\'Vk' > 

VdetiV- detiV+ kk' VCCt+^+|m| 

The above equation is not yet very useful because the states vk and uk also depend on the 
operator C. However, the overlap will not change if we perform unitary transformations (but with 
det/, det/^ = 1 rather than a phase; this restriction is meaningless because the single particle wave 
functions were defined only up to phase - more about this later) independently in the spaces Vl 
and Vr. This means that the bases of eigenvectors, {'Uk} and {wk}, can be replaced by arbitrary 
orthonormal bases of Vl and Vr. All the dependence on C becomes now explicit and we have 
formally established the connection 

detC(A)o A< i- |i+ >A (2.18) 

By an identical line of reasoning we can also arrive at the following relation connecting det C^ and 
the overlap of the many body states |i2± >a obtained by filling all the X^ states. 

A<R- \R+ >A = det(< Vk-|V'k^ >) 

KK' 

1 , 1 +, (2.19) 

= . = det < Vk , : C"\uk' > 

VdetN-detN+ kk' VCtC + m^ + \m\ 

Note that |i2± > represent the Dirac sea for — H^. We have the other formal connection 

detC^A)-^ A< R-\R+>A ■ (2.20) 

Comparing (2.17) and (2.19) it is clear that indeed (2.18) and (2.20) are conjugates of each other. 

The overlap is an adequate replacement for det C since, for states Vi for which the norm of 
Cvi is small relative to m, all the additional factors that appear in the overlap are constants. In a 
good choice of gauge fixing the smallness of the norm of Cvi means, for small gauge fields, that the 
typical momentum modes contained in the state Vi are small relative to m. We anyhow will have 
to regularize the expressions to make them well defined and any regularization induces mutilations 
of high momentum modes. Therefore, the fact that the overlap has not come out exactly det C 
even at the formal level is of no concern. 

Formally, the expression is also gauge invariant. Let Ql{G) denote the unitary representation 
of a gauge group element G in Vl and Qr{G) denote the representation of the same gauge group 
element G in Vr. The structure of C is such that one assumes [QR[G)yC[A^)QL{G) = C[A^); 
then the overlap is gauge invariant. The phases of det^j < li)^ (C/^(G))^ii)^ > and of det^j < 
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.,(^)n,fn-\Ji) 



v^j^' <3 l{G)v^£ > cancel in the expression below: 



det I < u^S ; n^^)^ > 

C(A«)(C(A«))t + m2 + |m| 

det ( < u^S , ^ ^^.KVl > I ^^'^^^ 



det < u^^(^«(G))tu^^ > det < ^^'^^^(GK^ > 

I'j \ J I'j \ 

The above equation holds because the sets {■u^! } S'^d {'u)^ } of the spaces Vj^ and F/^ respectively 
are complete. 

The overlap is still only a formal expression representing det C. We know that for small and 
slowly varying gauge potentials (with respect to the mass m) we would, in perturbation theory, 
get to leading order in — expressions identical to the possibly divergent, ordinary, unregulated 
Feynman diagrams. Thus we have not lost anything by replacing the previous formal expression 
for the chiral determinant by the overlap. We claim that we have made even some progress: It will 
be easier to turn the overlap into a well defined object rather than det C. Moreover, the formulation 
is such that we will not be tempted to consider an eigenvalue problem for C and as a result break 
Lorentz invariance. 

We can easily find gauge invariant ways to regulate the Hamiltonian problems. This may 
come as surprise because a gauge invariant regularization of a chiral anomalous theory must be 
pathological. When H gets regularized the simple structure it has in terms of the split V = Vr®Vl 
will be mangled. More precisely, with respect to the split of V we had: 

H±=(ct _B±) [B^C] = (2.22) 

Upon regularization, the commutativity made explicit above will be lost. We cannot therefore use 
the explicit construction of the overlap in terms of the two independent bases given earlier. As 
a result, in particular, the "proof" of gauge invariance of the overlap may break down. We still 
can easily define the second quantized system and say that we need the overlap between the two 
ground states. This defines the absolute value of the overlap but not its phase. The Hamiltonians 
can be regularized in a gauge invariant manner but the overlap is not guaranteed to come out gauge 
invariant. Physically it is clear that the control of "high momentum modes", the usual objective of 
ultraviolet regularization, is attained by regularizing the Hamiltonians. The possibility of anomalies 
is kept open by the existence of an additional ambiguity in the phase choices of the two second 
quantized ground states. The fact that usually one associates anomalies with the regularization is 
misleading. If the single ambiguity available to the theory is the fact that it has an uncontrolled 
ultraviolet behavior then anomalies are forced to use that opening to creep out. The overlap 
formulation is making it possible to deliberately leave an ambiguity in the regularized theory so 
that the anomalies can be reflected by it. In this way the problem of anomalies is separated from 
the problem of regularization and the truly ultraviolet infinities can be controlled by standard 
means. In short, the overlap expression provides a route to particular regularizations that break 
gauge invariance in a sort of minimal way. The obstacles to making these regularizations non- 
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perturbative (i.e. having them defined for any vector potentials A, without any restrictions, like, 
for example, on the topological sector) disappear. 



3. Zeros in nontrivial backgrounds 

If the gauge background cannot be deformed to zero, it has non-zero topological charge and 
the chiral determinant vanishes. Thus, for most physics questions, we can do without a definition 
of the phase of the two second quantized ground states of the Hamiltonians for gauge fields not 
deformable to zero. Before regularization, at the formal level, it is easy to see that the overlap 
will vanish when the operators C or C^^ have zero eigenstates. Consider a gauge field of instanton 
number one. C has a zero mode (Ci;o = 0) but C^^ has no zero mode. For all \\^ \ > m, the 
solutions vk and uk in (2.14) are paired. But there is one solution vq with \^ = — m that is 
not paired. This implies that there is one more ipj^ in (2.15) than there are V'k" iii (2.16) and 
the overlap matrix < 'ipj^~'ipj^J > is rectangular. The determinant of a rectangular matrix vanishes 
identically (the determinant is a function of the rows (columns) of a matrix which is linear in 
each and totally antisymmetric; if the matrix is rectangular there aren't enough columns (rows) 
for a non-vanishing totally antisymmetric function to exist). In a second quantized formalism the 
vanishing is attributed to the fact that one of the vacua has less filled states than the other. After 
regularization this argument will be more difficult to see because of the loss of the simple structure 
of H^. Nevertheless, the mechanism for such a vanishing to occur is sufficiently general that it 
can survive: All that has to happen is that the number of negative energy levels for one sign of 
the mass term be different from the number of negative energy levels for the other sign of the mass 
term. Even if the detailed simple connection between C and H^ is lost upon regularization, the 
departure from the unregularized structure is small for slowly varying states and the imbalance in 
the number of filled states should be preserved. To alter the imbalance some low momentum state 
would have to travel a large distance when the regularization is turned on because for any sign of 
the mass there is a large gap (> 2m) in the spectrum of H . 



4. The phase of the Chiral determinant 

Upon regularization, it is necessary to define the phase of the chiral determinant in the zero 
topological sector. This definition is arrived at by completely defining the two ground states making 
up the overlap. These states are parametrically dependent on the background vector potential. The 
phase of the overlap is the phase of a reasonably generic complex valued function defined over the 
space of vector potentials. There are points in that space where the overlap will vanish (these are 
"accidental" vanishing points - we are in the trivial topological sector). Generically, the phase 
will non-trivially wind around the vortex in the space of gauge potentials where the overlap stays 
zero. This is important because the inevitability of anomalies in the continuum in the class of 
acceptable regularizations can be established by showing that there are discs in the space of vector 
potentials whose boundaries live completely on a single gauge orbit, but, nevertheless, the phase 



of the chiral determinant must non-trivially wind when one goes around the boundary, since the 
chiral determinant has a zero of generic type in the interior of the disc. It would therefore be 
wrong to replace the definition of the chiral determinant with a definition that treats the absolute 
value and the phase completely independently. Upon regularization, such a definition might lead 
to an expression where the phases wind around points where the absolute value does not vanish, 
and points where the absolute value vanishes but the phase does not wind around. Therefore, the 
overlap, if properly defined as a complex valued function over the space of vector potentials, has 
the potential to preserve the topological understanding of anomalies in the continuum [10]. 

We wish to abstract from the formal overlap some properties about its phase that can survive 
regularization and the associated distortion of the structure of H^ relative to the split V = Vr®Vl. 
Assuming no accidental degeneracies the rays to which the states |X± >a and |i2± >a correspond 
are uniquely defined. To pick a state along the ray we need one additional condition. One way to 
single out a state is to fix the phase of its overlap with some fixed reference states I^Li > and 
|r2/^± >. The state \^l^ > needs to be "transverse" to the ray e'^^lXi >a, i-e. < f^L i l-^i >A7^ 
and then we can impose < Ql i l-^i >A> 0.* Similarly, we can impose < Qr ± |i2± >a> 0. It is 
important that the choice preserves the following properties of the formal determinant: 

1. The chiral determinants of opposite chiralities are conjugates of each other, i.e., 

detC(A)= [detC^A)]*. 

2. For a gauge potential A^^^x) and its Lorentz transformed partner A' Jx) = Ai^^A^[A~^x), the 
chiral determinants are formally related by 

det(C(A)) = det(C(A')). 

3. For a gauge potential A/^i^x) and the partner under parity 

^U^;) = {-A{-^>Xi)>^i(.-^>Xi)}> 

(in two dimensions X4 should be replaced by X2 above) the chiral determinants are formally 
related by 

[det(C(A'))]* = det(C(A)). 

4. For a gauge potential A/^i^x) and its partner under charge conjugation A' Jx) = —A*Jx), the 
chiral determinants are formally related in 2-D by 

det(C(A')) = det(C(A)), 

and in 4-D by** 

[det(C(A'))]* = det(C(A)). 



* It is possible that the overlap with the reference state vanishes accidentally and our definition 
for the phase breaks down. We believe that the gauge potentials where these accidents happen can 
be safely ignored because the reference states are not that special. 

** This property has an important consequence in 4-D. The chiral determinant will be real when 
A is restricted to a real or pseudoreal representation (in which case there exists a matrix in the 
internal space R such that A' = RAR~^). 
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5. For a gauge potential ^^(cc) and its partner under a global gauge transformation A' {x) = 
g^A^[x)g where g G U(N), the chiral determinants are formally related by 

det(C(A')) = det(C(A)). 

We show below in a sequence of Lemmas that the Wigner-Brillouin phase choice for the many 
body states preserves the above properties of the chiral determinant when written as an overlap 
of these two many body states. The Wigner-Brillouin phase choice is simply demanding that the 
overlaps o< -K ± |-ff± >a and o< L i: |X± >a are real and positive. The important point for us 
is that only quite general properties of H^ are used in establishing 1-5 for the overlap. Even if 
regularization spoils the simple structure of H^ (which has to happen as we explained before), the 
properties of the overlap under the symmetries can be preserved because the Hamiltonians will still 
obey the general properties used in proving the Lemmas. Of course, for a lattice regularization one 
must give up on Lorentz invariance, but one still has the usual invariance under the appropriate 
discrete crystallographic group. 

It is worthwhile to insist on the discrete symmetries: For example, in four dimensions, the 
combined requirements of translational, Lorentz and global gauge covariance together with the 
above behavior under parity and charge conjugation imply that the phase of the overlap, when 
expanded in powers of the vector potential A^ starts at third power. This is a desirable property 
for a regularization to have. 

Now we proceed to show that 1-5 hold for the overlap formula with the Wigner-Brillouin phase 
choice for the many body states. Essentially, 1-5 are consequences of the fact that |i2± >o and 
|£± >o are invariant under the symmetries under consideration. 

In what follows the following Lemma about a general unitary matrix will be useful. 

Lemma 4.0 Let 

\C D 

be an N X N unitary matrix with A being a K X K square matrix (0 < K < N). Then 

detA 



detU 



detDt 



Proof: 

We write 

det U = detA det{D - CA'^B). 



For a unitary matrix, one has the relation 

D-CA-'B = D 



iR - nt-i 



Inserting this in the previous equation proves the Lemma ■ 

The single particle eigenfunctions for the gauge free case are chosen to satisfy the following 
two properties: 
1. Since 

-['y4]a^il^(.(.-x,X4)'yi,{-y,y4)Sj;0)['y4]sf3 = il^{xai,yPj;0) (4.1) 
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we set 

'iP^^ixai; 0) = [liUi^^^ii-x, x,)pi; 0) (4.2) 

74 is replaced by 72 in 2-D. 
2. Given the single particle states for the A = problem let U{+m; 0| — m; 0) be the unitary matrix 
connecting the basis at +m to the basis at — m. We assume that the single particle states have 
been chosen so that o< R — \R+ >o is real and positive. Then from (4.2), o< L — \L+ >o is 
also real and positive. From Lemma 4.0, we then have 

detZ^(-m;0| +m;0) = 1 (4.3) 

Lemma 4.1 If o< -K ± |-ff± >a and o< L i: |X± >a are reai and positive then 

A< R- \R+ >A= A< L+ \L- >A 



Proof: 

From Lemma 4.0, we have 

detZ^(±m;0|±m;A) = 44^^|^. (4.4) 

A< L ± \L± >o 

Since o< -K ± |-ff± >a and o< L i: |X± >a are real and positive (4.4) implies that 

detZ^(±m;0|±m;A) = 1 (4.5) 

Since 

U{-m; A\ +m;A) = U{-m; A\ - m; 0)U{-m; 0| + m; 0)U{+m; 0| + m; A), (4.6) 

(4.3) and (4.5) implies that 

det U{-m; A\ + m;A)=l (4.7) 

By Lemma 4.0, we have 

det U{-m; A\ + m; A) = a<R-\R+>a ^^_g^ 

A< L + \L- >A 

and this along with (4.7) proves the Lemma ■ 

Lemma 4.2 For a gauge potential A/^i^x) let its Lorentz transformed partner be 

A'^{x) = A^,A,{A-'x). 

If o< -ff i l-ffi >A and o< -K ± |-ff± >a' are reai and positive then 

A'< R — \R+ >A'= A< R — I-K+ >A ■ 
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Proof: 

Since, 



Yi.^[xai,yPJ-A'] 



Ju>(A) 



Ci^ 



^^[{K-^x)^i,{K-^y)6j-A] 



-iu>(K) 



5/3 



(4.9) 



g8w(A)^^g 8w(A)^ ^g have the following relation between the eigenvectors of the two 



where 7^A^i. 

Hamiltonians (assuming no degeneracies): 



iP'^^{xai-A') = e 



^4>T ('e'^(^) 



afi 



i^^^{K-^xPi-A) 



(4.10) 



In the above equation ^j^ are undetermined phases. For the case when A = 0, (4.9) implies 
degeneracies and therefore 



iP'^^{xai-Q)= e 



,«^(A) 



a/3 



iP^^{K-^xf3i-Q) 



(4.11) 



defines another set of eigenvectors for the free Hamiltonians. The transformation matrix relating 
the vectors tf}^ and t^j^ will be block diagonal and will only mix vectors in a degenerate subspace. 
Further, the individual blocks will all have unit determinant. Since the state |i2± > is obtained 
by filling all states in a given degenerate subspace, using tf}^ or t^j^ will result in the same state 
\R± >. From (4.10) and (4.11) we have 

o<R± \R± >A' =o<R± \R± >A e Sa- *a-* (4.12) 

We are given that o< -K ± |-ff± >a and o< -K ± |-ff± >a' are real and positive and therefore we have 



K 







(4.13) 



From (4.10) we have 



=»E..(^f-"-^r) 



(4.14) 



A'< R - \R+ >A'= A< R- \R+ >A e 

and this along with (4.13) proves the Lemma ■ 

Lemma 4.3 For a gauge potential A/^i^x) and the partner under parity 

^U^;) = {-A{-x,X4),A4{-x,X4)} 

(in two dimensions X4 should be replaced by X2 above), if < -K i l-Ki >a, < -^ i l-^i >A, 
o< -ff ± |-ff± >A' 3.nd o< i ± |X± >A' are reai and positive then 



,< R- \R+ >, 



i< L- \L+ >A . 



Proof: 

For the gauge potential A/^i^x) and its partner under parity. A' Jx) we have the relation 

ll'^[xai,y(3j;A'] = -[j4]aj'H.'^[{-x,X4)ji,{-x,X4)6j; A][j4\si3 
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(4.15) 



In the above equation, 74 should be replaced by 72 in 2-D and the proof that follows trivially 
extends to 2-D. The eigenstates for A' are related to the eigenstates for A by 



,^K- 



R±, 



iP'^^{xai-A') = e'^K [74]„^Vk (-^a;4/3i; A) 



^l^ixai; A!) = e'^K [j^]^^i;^^{-xx^(3i; A) 
In the above equations, 4>j^ and 4>j^ are undetermined phases. Using (4.2) and (4.16) we get 

o<L± \L± >A' =o<R± \R± >A e ^K ^i^ 



(4.16) 



o<R± \R± >A' =o<L± \L± >A e ^K *A-* 



(4.17) 



Since o< -ff ± \R± >a, o< L± \L± >a, o< -K± |-ff± >a' and o< L± \L± >a' are real and positive, 
it follows that 

E^x^ = E^x^=0 (4-18) 



From (4.16) we have 



!-< R - \R+ >A'= A<L- \L+ >A e^K^'^^^-'^K ) 



(4.19) 



and this along with (4.18) proves the Lemma ■ 

Lemma 4.4 For a gauge potential A/^i^x) and its partner under charge conjugation A' Jx) = 

-A*^{x), if o< R ± |-ff± >A, o< L ± \L± >A, o< -ff ± \R± >A' and o< L ± \L± >a' are reai and 
positive then 

A'< R — I-K+ >A'= A< R — I-K+ >A 

in 2-D and 

A'< R - \R+ >A'= A< L- \L+ >A 

in 4-D. 

Proof: 

For the gauge potential A/^i^x) and its partner under charge conjugation, A'Jx), we have the 

following relation in 2-D 



Il^[xai,yPj;A'] = - ['yi]a^'ii^[x'yi,ySj;A]['yi]sp 
and the following relation in 4-D 



ll^[xai,y(3j;A'] = [S2]„^H±[a;7i,2/,^i; A][S 



2J5/3 



where 



(72 

-(72 

In 2-D, (4.20) gives (assuming no degeneracies) 



i^^H^^az; A') = e<^ [j^U [i'K^] *(^/3^; ^) 
,l;^,^{xaz;A') = e<^ [j,U['^^,^]\xPz; A) 



(4.20) 

(4.21) 
(4.22) 

(4.23) 
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In the above equations, 4>^ and 4>^ are undetermined phases. For the case when A = 0, the 
relations corresponding to (4.23) are 

'iPf^ixaz; 0) = e^^K [^,]„^ [^^±] \a^p^. Q) 
'iP^^ixai; 0) = e^^K [^,]„^ [^£.±] *(a;/3i; 0) 

In the above equations, Xk ^^^ undetermined phases. From (4.24) we have 

o<R-\R+ >o= [o<L- |X+>o]*e'SA-(>^A--XK) (■4_25) 

The second property of the many body states in the free case, namely o < R — \R+ >o and 
o< L — \L+ >o are real and positive, implies that 

Y.'^xti - Xji) = (4.26) 

K 

From (4.23) and (4.24) we have 



(4.27) 



Since o< -ff ± \R± >a' , o< L± \L± >a' , o< -K ± |-ff± >a and o< L± \L± >a are real and positive, 
it follows from (4.27) that 

E^"" = E^a" = ExI (4.28) 

K K K 

From (4.23) we have 

A'< R - \R+ >A' =[a<L- \L+ >^]*eEK(*A-+-*A--) 



A<R-\R+>Ae^K(^^^^-'^^^~^ 

A< R- \R+ >A C^k^^^-^k) 

A< R — \R+ >A 



(4.29) 



which proves the Lemma in 2-D. In the above equation we have used Lemma 4.0 in the second 
equality, (4.28) in the third and (4.26) in the last equality. 
In 4-D (4.21) gives 

',p^^{xai;A') = e^^K^[i:2U['>P^^]\xPi;A) (4.30) 

where 4>j^ are arbitrary phases. For the case when A = 0, we have the situation where 

'>P^^{xai;0) = e'>^"'[S2]„4Vr]*(a;/3i;0) (4.31) 

are another set of eigenvectors. Xk ^^^ phases chosen in such a way that the many body states 
|i2± > formed out oiip^ and ip^ are the same. From (4.31) we have the relation 

o<R- \R+ >o= [o< R - \R+ >o]*e'SA-(>^A-+->^A-") (4.32) 
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We have chosen the states |i2± >o so that o< R — \R+ >o is real (c.f. second property of the free 
states) and therefore it follows that 

Exr = Exr (4-33) 

K K 

From (4.30) and (4.31) we have 

o< ii;± \R± >A'=o< R± \R± >A*e'SA-(*A-*->^A-*) (4.34) 

and since o< -K ± |-ff± >a' and o< -K ± |-ff± >a are real, we get 

E^"" = Ex^^" (4-35) 

K K 

From (4.30) we also have 



v< R - \R+ >A' =[a<R- \R+ >A]*e^K(^K^-^K~ 

E/ R+ _ R—\ 
= A<L- |X+ >A 



) 



(4.36) 



where we have used Lemma 4.0 in the second equality, (4.35) in the third and (4.33) in the last 
equality. This proves the Lemma in 4-D ■ 

Lemma 4.5 For a gauge potential A/^i^x) and its partner under a global gauge transformation 
A' Jx) = g^Ai^[x)g where g G U(N), if o< R ± |-ff± >a and o< -ff ± |-ff± >a' a.re reai and positive 
then 

A'< R — I-K+ >A'= A< R — I-K+ >A 

Proof: 

For the gauge potential A/^i^x) and its partner under a global gauge transformation. A' Jx) we 
have 

ll^{xai,y(3j;A') = [g%ll^{xak,y(3l; A)[g]i, (4.37) 

As a consequence of this we have the following relation 

i;^,^{xai;A') = e^'^"^ [g^],,i;^,^{xaj; A) (4.38) 

For A = we have degeneracies and therefore 

i^^^ixai; 0) = [g^],,i^§^{xaj; 0) (4.39) 

is another set of eigenvectors for the free Hamiltonians. The transformation matrix relating ip^ 
and tpj^ will be block diagonal and will mix only states in a degenerate subspace. The states 
|-ff± >o are obtained by filling all the states in a subspace. From (4.38) and (4.39) we have 

o<R±\R± >A'=o<R±\R± >A e'SA-(*A-*+«) (4.40) 
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where e'" = [det^fJN Since o< -K i l-Ki >a and o< -K ± |-ff± >a' are real and positive it follows 
from (4.40) that 

E^"" = E" (4-41) 



From (4.38) and (4.41) we have 



!-< R - \R+ >A' = A<R- \R+ >A eSA-(*A-+-*A- ) 



'A' — A^ -f- - l-a-l- ^A 

A< R — \R+ >A 



(4.42) 



and this proves the Lemma 



5. Expectation values of fermionic operators 

In section 2, we obtained a formula for the chiral determinant as an overlap of two many body 
states. To the chiral Dirac operator C(A) in the presence of an arbitrary hermitian matrix valued 
vector potential A, two single particle Hamiltonians H^(A) given by 

were associated. 

H^{A)= Y, al^^Il^{xai,yPJ;A)ayp, (5.2) 

xai,yf3j 

are two many body Hamiltonians for non-interacting fermions with the single particle Hamiltonians 

'''x/di 



given by (5.1). a^ai and a^^^ are single particle destruction and creation operators obeying canonical 



commutation relations 

{axai,ayi3j} = 0; {a:,ai,ayf3j} = 0; {a-xai^a-ypj] = ^xy^aP^ij- (5.3) 

It was shown that 



det C{A) <:^ A< L - |X+ >A 
detC^(A)o A< R- \R+ >A 



(5.4) 



were proper definitions of the chiral determinants. The |X± >a were ground states oiTi {A) and 
the |i2± >A were ground states of —7i {A). Explicitly, referring to (2.11) which defines the single 
particle eigenfunctions of H (A), if we define new creation operators 



K^] = E «L.V'r (^"^; A); [bf^] ' = E «L.V'^'^(^"^; a), (5.5) 



then 

Li?±lt|n^. ir_L ^_ TTrAl^ilti 



i^± >= nt^K^] > >■' i^± >= nt^K^] > > (5-6) 



K K 

16 



where |0 > is the bare vacuum annihilated by all the destruction operators a^ai- If we define 

xai 
xai 

OkUA) = Y,[i^k-i^ocz;A)Yi;^t{xaz;A) 

xai 



(5.7) 



then the overlaps are given by 



LL/ 



< R-\R+ >A= det O^^(A); a< L - \L+ >a= det O^^(A). 



(5.8) 



We wish to construct generating functionals for the expectation values of physical fermionic 
operators in an arbitrary gauge background. If the gauge background is trivial topologically the 
generating functional will factorize into a product of generating functionals, one factor for each 
gauge group multiplet. Then one can construct each of these factors independently by replacing 
the formal path integrals 






*«D^ + ^^ 



*^D^-^^ 



c\ ' / j Kixai xai ' Ixai xaiJ 

xai 

c\ ' / j \ Ixai xai I Ixai xaij 



by well defined objects, Z^[fj^,r]^; A) and Z^[fj^,r]^; A) respectively. ^^ and ^^ are four com- 
ponent Dirac spinors that represent the physical two component right handed and left handed 
fermions and the terms bilinear in them are the two chiral components of the Dirac operator as 
defined in (2.3). The bilinear terms contain chiral projectors ( — j^) and are therefore not invert- 
ible. This implies that in the first (second) of the two equations above only half the modes are 
propagating and represent a physical single right (left) handed fermion. rj^, rj^, fj^ and fj^ are 
Grassmann variables and are introduced in the usual manner to facilitate the definition of fermionic 
expectation values. The expressions for Z^ and Z^ are 



Z^{fi^, r]^; A) = A<R-\ exp [X](7?^„,a,„. + VxaiaLi) 

xai 

Z^iv^, V^; A) = A<L-\exp\ ^{fj^c^.a^ai + 7?L.«Lj 



\R+>A 
\L+>A 



(5.9) 



Completely antisymmetrized re-point functions of fermions are obtained by differentiation with re- 
spect to the associated Grassmann variables. Note that the dependence on fj^ ~^^ , fj^ ^ ^ ; V^ 2 ^ 
and r]^ ~^^ is no longer as trivial as it was in the formal path integral expressions. Nevertheless, 
we have not added any extra physical modes: Consider an re-point function and order the operators 
inside so that all the annihilation operators appear to the left and the creation operators appear 



17 



to the right. In doing this there will be many terms with either n or less operator insertions. Ex- 



6-± 



b^^ 



>^K^>^K^ basis. Then, 



propagate. Although we are 



press the fermion creation and annihilation operators in the 

from (5.5) and (5.9) it becomes clear that only half of the modes wil 

representing a single right (left) handed fermion by a four component Dirac spinor we do have the 

correct number of physical propagating modes. The explicit formulae in subsection (5.1) will make 

this clear. 

If the gauge background is not trivial topologically the fermionic functional integral of the 
source term no longer factorizes and we must look for a definition that includes all the physical 
fermionic fields in the theory together. Eq. (5.9) generalizes easily, multiplying the expressions 
for each gauge multiplet. The product of the vacua matrix elements can be written as the matrix 
element of a complete source term between the two full vacuum states, each given by a tensor 
product over the multiplet vacua: 



^< 0- |exp[ Y. (^ri"af^^« +7?,^i«(a^'^«)L)+ E iv^a^t'f + V^a'na^''nl^.)\\0+ >A 

|0± >A = (g) \R, Ir; ± >a (g) \L, II; ± >A 

Ir II 

(5.10) 
The remaining task is to define the phases of the two complete vacua |0± >a- This will be done 
in section (5.2). 

5.1. Zero topological sector 

In the zero topological sector one can write Z^ and Z^ in (5.9) as 



Z^(7?^, 7?^; A) = exp[Ty^(7?^, t?^; A)U< R - \R+ >a 
Z'^ifj'^, rj"-; A) = exp[Ty^(7?^, 7?^; A)U< L - \L+ >a 



(5.11) 



The states |i2± >a and |X± >a are well defined since their phases have been fixed by the Wigner- 
Brillouin choice (that is demanding that their overlap with the corresponding free states be real 
and positive). In particular, the two point functions are fully defined and are 



< '^xai'^yPJ >A O 



< '^xai'^yPJ >A O 



A< R — \R-\- >A 
t 

xat 



S/3j4«»)I-^+ >A 



(5.12) 



< L-\L+ >, 



By referring to the formal path integrals we see that these are the inverses of the chiral operators 
C^ and C respectively. Therefore, we would expect to have 



< ^xai^yPj >A- 



^ ^yPJ^xai >^ 



(5.13) 



We proceed to show that the precisely defined right hand sides of (5.12) also satisfy (5.13). We use 
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(5.5) to write a^ai and aj,^^ as 



'i" 



(5.14) 



From (5.6) it follows that 



[b'^+y\R+>A=0; b'^+\R+>A=0; A<R-K~ = 0; A< R - {[b'^-y = 0; (5.15) 



Then 



-A< R - \R+ >A ^xy^aP^ij " -A< R " K^LiS/^i " °'yPJ°'lai)\R+ >A 



A 



A 



--A< R - \ay[3ja^^^\R+ >A 



E V'ii(2//3j)[V'^^r(^az)j <Q\x{b^K-b^,-y^\\ npr 



Ki ,K'^ K' 



A< L - \ayp-axai\L-\- >a 



A 



"i^IlH^ 



10 > 



|0> 



-A< -Z^ - l(%/3,-^a;ai " axaiay0j)\L+ >A +7;A< L - \L+ >a ^xy^ajd^i. 



(5.16) 



We have used (5.3) in the first and last equality and (5.14) and (5.15) in the second and fourth 
equality. The matrix elements appearing after the second equality are overlaps bewteen two states, 
each obtained by filling all the positive energy states and one negative energy state for both Hamil- 
tonians respectively. The matrix elements appearing after the third equality are overlaps between 
two states, each obtained by filling all the negative energy states but one for each of the two Hamil- 
tonians. There is a one to one correspondence between the two overlaps given by a particle hole 
transformation. More specifically, the two overlaps are related by Lemma (4.0): 



detZ^(-m; A| +m;A) 



<0|nx, 6^.76^.7 



^i^t 



n, 



6-+ 



10 > 



< 0| Ha- &x" 



A 



&i^t n, 



h^<^ 



10 > 



(5.17) 



By (4.7) the left hand side above is unity and this establishes the third equality in (5.16). Dividing 
(5.16) by A< R — \R-\- >A= {a< L — \L-\- >a ) we obtain the analogue of (5.13). 

An explicit expression for the two-point function can be obtained by writing the determinants 
representing the matrix elements in (5.16) in terms of the matrices O defined in (5.7). We employ 
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the identity det I ] = detXdet(F — ZX ^Y) where now y is a column vector and Z is a 

row vector and F is a number in the second equality in (5.16) to obtain: 



A< R — \R-\- >A 



Ki ,K[ 



qLL _ 'Sp qLR 



o 



RR 



K,K' 



pRL 
KK' '''''- 



(5.18) 

Using the explicit expressions in (5.7) and the completeness relation for the eigenvectors of H^(A), 
(5.18) reduces to 



A< R — \R-\- >A 



Y, i^^tivPj; A)[0n7':j, \i^§-{xaz; A) 



(5.19) 



KK' 



The right side of the above formula has a simple interpretation. Since we are computing the 
propagator of a right handed fermion we expect to get [0^^]~^ . The two wavefunctions carry 
out a change of bases from the two bases of the negative energy subspaces of H^ to a common 
coordinate basis. It is evident that the matrix O^^ is treated as a representative of a mapping 
between two distinct spaces. The right hand side viewed as a matrix element implies that the 
associated matrix has a rank equal to the dimension of the space spanned by the index K which is 
nothing but the basis for a single right handed fermion. 

(5.19) can be recast in a form that makes the relation (5.13) transparent: Since 



/qRR qRL\ 



is a unitary matrix it follows that 



O 



RR 



o 



RR 



o 



LR 



o 



LL 



o 



RL 



(5.20) 



Employing again completeness we obtain the following identity: 



E 

KK' 



i;^t{yPJ;A)[0^^]-],, i;^-{xai;A) +i;^-{yPJ;A) {O^^Y i^f^{xai;A) 

L J L J KK' L 



xy^af3^ij 



(5.21) 



Using (5.21) the expression for the two-point function in (5.19) can be rewritten as 
1a<R- \{aL^ayP3 - S/3j«LJI^+ >a 



i<R- \R+ >, 



E 

KK' L 



il;^+{yPJ; A)[0''%,],, Vx" (^"^i ^) " V'F (j/Z^i; A) (O^^)^ il;f+{xai; A) 

KK' 



(5.22) 

To obtain the two-point function for the left-handed fermion, all the R's have to be replaced by L's 
in the above equation and the relation (5.13) is made transparent. 
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We will presently show that W^^fj^jT]^; A) is bilinear and has only terms connecting fj^ and 
7/^. The same is true for the left handed fermions. This has the following consequences: The 
TC-point functions should have an equal number of fermion creation and annihilation operators for 
them to be non-zero. The re-point functions are sums of products of the two-point function in (5.12). 
The proof is straightforward but slightly complicated because the expectation value is between two 
different states |i2± >. From (5.5), we use 



K-'y= Y.^L.i^KH^<^^■,Ay, [bffY = ^aLV'^.-(^az; A), 



(5.23) 



fi+it 



uL-lt 



to express all the a].^- in terms of [bj^ ] and [bj^ ] and we use 

r -| * r 

bff = "^a^ai V'K"(a;ai;A) ; 6^.+ = ^ a^„i tpj^'^ixai; A) 



(5.24) 



to express all the a^ai in terms of bj^~ and bj^~ . The above transformations are non-singular (and 
therefore valid) if the overlap a< R—\R+ >a is non-zero. One way to see this is as follows. In (5.23) 
we know that all the 'ipj^'^[xai; A) are linearly independent and so are all the 'ipj^~[xai; A). All we 
need to show is that we cannot write any 'ipj^'^[xai; A) as a linear combination of 'ipj^~[xai; A). If 
this were possible then b^~\R— >a= and we would have a< R— I-K+ >a= 0. Similar arguments 
can be used for (5.24) also. Therefore we can write 



a^ca^ = Y, ^K (2;"^; ^)^X + Yl <^K^(=^"^; ^)^ 



K 



aL, = Y\cl>^-'{xaz;A) [b^+Y + ^Uff {xaz; A) [bff] 



LL-lt 



(5.25) 



where 4>^ [xai; A) and 4>^ [xai; A) are obtained by inverting the relations (5.23) and (5.24). That 



IS 



^</.^.-(a;ai; A)[V'x-(2//3i; A)j* + ^</.^+(a;ai; A)[V'A+(2//3i;^) 



^xy^a(3^ij 



Y\ci>f^{xai-A) Vx+(2//3i;A) + ^U^r(a;ai;A) i^^f{yf3j-A) = 6,,y6^,p6 



(5.26) 



Let 



Q- 



Sh 



E 

xai;K 

E 



fj^^,cl>^-{xai; A)b^- + 7?^„, cl>^f{xai; A) 



\L+ 



fi^,,(t^^^{xai;A)b^^ + r]^,, ct^^-' {xai; A) 



R aR+i 



b^- 



(.«+lt 



xa.i;K _ 

From (5.3), (5.23), (5.24) and (5.26) it follows that the commutator 



(5.27) 



[Q+,Q-]= Y vL^VyPjY. 
xai;yf3j K 



'R 



</.«- {xai) V'r {yPJ) + i;^-^{xai) ct>'^+ {y(5j) 



R+i 



kR+i 



7 y VxaiVyPj 7 ^ 
xai;yf3j K 



\L+ 



cl)f+{xai) VF(2//3i) +il;^-{xai) cl)ff{yPJ) 



^xy^a(3^ij 
^xy^a(3^ij 

(5.28) 
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is a c-number. There are only terms that mix rj^ and rj^ in the commutator. All other terms are 
zero because a^^j = and [aj.„j]^ = 0. Upon inserting (5.25) in (5.9) and using (5.27) we see that 

exp[W^{fj^,rj^; A)]a< R - \R+ >A = A< R - | exp[Q_ + Q+]|ii;+ >a 

= exp(^[S+, S-])a< R - I exp[S_] exp[S+]|ii;+ >a 

= exp(-[S+, Q-])a< R - \R+ >A 

(5.29) 
We have split the exponential in the second equality and used the fact that the commutator in 
(5.28) is a c-number. In the third equality we have used (5.15). Therefore, 

Ty^(7?^7?^;A) = ^[S+,S_] (5.30) 

where the expression for the right hand side is in (5.28) and is a bilinear in t/^ and t)^. 

5.2. Non-zero topological sectors 

Our construction for nonzero topological charge will apply only to situations where a local fermionic 
operator J^ exists with the following properties: It is a Lorentz scalar; it is translational and gauge 
invariant; it can connect a state |0+ >a to a state |0— >a when A carries topological charge. 
We shall take J^ to be the most local and least dimension operator satisfying these requirements. 
Such an operator will be referred to as a 't Hooft [11] vertex since it would exhibit the (minimal) 
violation of global symmetries due to instanton-like gauge fields. If such an operator does not 
exist we have no physical reason to try to define fermionic expectation values in nontrivial gauge 
backgrounds. Typically, there will be one or few such operators. They make up a fundamental 
set in the sense that any nontrivial gauge background that would give a nonvanishing expectation 
value to some fermionic operator will also give a nonvanishing expectation value to a monomial 
made out of the fundamental 't Hooft vertices J^ and their hermitian conjugates. Henceforth we 
shall assume, mainly for notational simplicity, that there is only one such operator J^. 

We shall define the phases of the states |0± >a in (5.9) by replacing the physical fields in J^ by 
the matching creation and annihilation operators as in (5.9), and denoting the resulting operator 
by f demanding a< ± |J^"|0± >a> or a< ± |[J^t]"|0± >a> for some n > 0. One may 
unify our phase definitions to one equation that is satisfactory for all topological sectors: 



A<0± 



e^ + e^' 



|0± >o is real and positive. (5.31) 

In a topological sector with topological charge q, Z in (5.9) can be written as 

Z{fj, rj; A) = exp[Wq{fj, rj; A)]z{fj, rj; A). (5.32) 

Wq is bilinear and only terms connecting an rj to an fj appear, z is a monomial in a subset of the 
Grassmann sources fj and ?/ of a power increasing with \q\. z reflects the insertions needed to "soak 
up" the zero modes in path integral formulations. This structure can be derived by following the 
steps described in the previous sub-section and therefore we omit the details. 



22 



6. The chiral determinant in perturbation theory 

For a gauge background that can be smoothly deformed to zero we can use perturbation 
theory. Since it is generally easier to carry out perturbative computations using Feynman diagram 
techniques [3], rather than Schrodinger time independent perturbation theory (which could have 
been directly applied to the overlap [4]), we take a detour and first derive an Euclidean path integral 
expression for the overlap formula. Once we have a path integral expression, the Feynman rules 
can be easily written down. 

To write the overlap a< L — \L-\- >a as a path integral we proceed as follows. For generic 
perturbative gauge fields, we expect q< L ^ |X± >a to be non-zero. Since |X± >a are ground 
states of 7i {A) given by (5.2), we can write them as 

\L±>A= lim ' ^° (6.1) 

|£± >A so constructed satisfies the Wigner-Brillouin phase choice. From (6.1), the overlap can be 
written as 

A< L — i+ >A= lim — , = — , = (6.2) 

To convert the above expression to a path integral, we first write down the one particle 
Schrodinger wave equation with the Hamiltonian H (c.f. (2.10): 

i5tV = 75[±m + 7^(5^ + M^(a;))]V' (6.3) 

We view this as a five dimensional Dirac equation and use the five- dimensional Minkowski metric 
( + ,—,—,—,—). The five dimensional 7 matrices are taken as Fq = 75, F^ = —27^, /x = 1,2,3,4. 
We now want to recast the above equation in the form of the five dimensional Dirac equation. To 
this end we multiply by Fq and obtain: 

(iF»5„ - FM„ T m)V' = (6.4) 

Here a = 0, 1,2,3,4, Aq = and A^i^x) is independent of t. The Minkowski path integral corre- 
sponding to the second quantized system equivalent to the above Dirac equation and incorporating 
Fermi statistics is: 

[#][#]e"'-^'^'^"^'^""'^"^^"''' (6.5) 



We now go to five dimensional Euclidean space, replacing t by —is to obtain the Euclidean path 
integral 

[di}][dil}]e- -I ^^^'^'^^^^^^^'^^'^^'"'^'^ (6.6) 



Using (6.6) the three factors in (6.2) can be written as a path integral and the result is 

/'S=T\ - 

s = — T_ 



A< L - \L+ >A= lim 






(6.7) 
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where the s— dependent Euclidean Lagrangian is 

LE{i^,il^;A,m,s)= j[dx\'ij){x,s) ^r,ds + ^i^{di^ + iAi^{x)) + me{s) tl}{x,s). (6.8) 

Note that the Lagrangians in the denominator are in fact independent of s. This is because in one 
factor s is always negative and in the other factor s is always positive. This reflects the fact that 
one Hamitonian is 7i~ and the other is Ti'^ . On the other hand in the numerator the Lagrangian 
changes as s goes from the negative side to positive side and this accounts for the two Hamiltonians 
7i^ present in the numerator. 

(6.7) is still unregularized. For finite T± gauge invariance is broken by the boundaries, no 
matter what boundary conditions we pick. When the T± are taken to infinity this breaking dis- 
appears from the absolute magnitude of the expression; at the same time the free propagator is 
specified by requiring decay at infinite s. The phase of the expression is determined and will typ- 
ically come out gauge dependent. Even before considering regularization we need to understand 
how in perturbation theory the infinite s ranges do not induce additional infinities.* 

In perturbation theory one computes the logarithm of the overlap. At a fixed order one 
has three diagrams to consider: one from the numerator and two from the two factors in the 
denominator. Each of these diagrams has n vertices and each vertex has, among other labels, an s 

label, si, S2, s„. The integrals over the si variables are to be done in the intervals prescribed and 

the limits on T± of the three diagrams are taken in a correlated way. To obtain the rules directly 
in the large T± limit, change variables in each of the diagrams from si,S2,....s„ to a center of 
mass coordinate s and n—1 orthogonal linear combinations; the complicated bounds of integration 
for the latter simplify in the limit T± -^ oo and there the ranges of the relative variables are 
independent and run over the whole real axis. The integrals over the relative coordinates are to be 
performed first; as long as the center of mass is kept finite no divergences will appear from these 
integrals. These integrals can be done independently for each of the three diagrams. One is now 
left with three s integrals: one over s coming from the numerators and two over s± coming from the 
denominators. The square roots in the denominator attach a factor j to the s± integrals. However, 
the diagrams coming from the denominator, once the integrals over the relative coordinates have 
been carried out over their infinite ranges have no s± dependence left as the systems corresponding 
to the denominator are homogeneous. (This will be only approximately true if we smooth out the 
behavior of the mass term but this has no effect on the finiteness.) The factors of one half simply 
cancel against the factors of 2 in the ranges and the constant terms are subtracted, respectively, 
from the positive s and negative s portions of the s integrands rendering the corresponding integrals 
finite. Since the denominators are real the subtractions only affect the real part of the diagrams. 
The imaginary part of the diagrams generated by the numerator can be summed over s without 
encountering divergences, i.e. are well defined. 

The above defines the unregulated Feynman diagrams of the theory. One may ask whether 
the path integral representations, or, equivalently, the limiting expressions leading to them can be 
taken over to apply also for A fields that cannot be deformed to zero. The answer is emphatically 
no: The path integrals, the limiting expressions and the related Feynman rules are to be used only 
for perturbative A's. For example, if the background is an instanton, the limiting expressions would 
be saturated by an excited state on one of the sides rather than vanish as they should. So, only 



For a detailed discussion see recent work by S. Aoki and R. Levien [12]. 
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the overlap formula is to be taken outside perturbation theory. 



7. Regularizing the Overlap 

We turn now to the problem of regularizing the overlap formulae (2.18) and (2.20) on the 
lattice. We choose the lattice because on it the dynamics of the gauge fields can be defined non- 
perturbatively and we wish to extend this definition to chiral fermions. 

More precisely, our task is the following: Assume, for definiteness, that Euclidean spacetime is a 
torus. On it we define a gauge connection that may be in a nontrivial topological sector. We embed 
a finite hypercubic lattice in this torus and define link variables U^i^x) as the parallel transporters 
along the link connecting x io x-\- fi. We wish to define two Hamiltonian matrices H with indices 
as in (5.2) (except x and y now label sites and have a finite range) such that, following through 
with the construction of the overlap will produce a functional of the link variables which will vanish 
if the continuum connection was in a nontrivial sector and will agree with continuum perturbation 
theory in the topologically trivial sector. By this we mean that if the external connection is small 
and slowly varying on the lattice scale (i.e. the lattice is fine enough), a double expansion in it 
and its derivatives has the form one might have obtained by regularizing the determinant of the 
chiral Dirac operator by some acceptable continuum method. If the continuum connection is in 
a topologically nontrivial sector, the lattice overlap should vanish due to unequal fillings of the 
two Dirac seas, as explained in section 3. In addition we require that if we follow through with 
the construction of the generating functional of fermion field expectation values (eq. (5.9)) for 
the lattice Hamiltonians in an any topological sector, the correct continuum behavior is obtained 
when the external momenta of the fermi fields are small relative to the inverse lattice spacing. 
We also require to preserve the symmetries mentioned in section 4, more specifically, eq. (4.9) for 
the hypercubic space symmetry group, eq. (4.15), eq. (4.20) in 2-D, eq. (4.21) in 4-D and eq. 
(4.37). The task is to find local matrices H^ (in the sense that matrix elements connecting sites 
that are too far away vanish exponentially) satisfying the above requirements. These matrices are 
parametrically dependent on the link variables with no restrictions on them and so is the rest of 
the construction. Thus we obtain a non-perturbative definition. 

The free Hamiltonian matrices must be traceless by (4.15) and so must be the Hamiltoni- 
ans in the presence of parity invariant gauge potentials. It causes no complications to make the 
Hamiltonian matrices traceless for all gauge fields. Then the structure is as in eq. (2.22), 

however we shall see that we need to allow for [B , C] ^ 0. 

On the lattice C becomes a finite matrix of fixed dimension and C^ its adjoint under the 
ordinary inner product. Thus C and C^ have the same rank. Under a local gauge transformation 
C will transform by conjugation and therefore det C would be gauge invariant. These properties 
show that det C would be an incorrect replacement for the chiral determinant since there would be 
no anomalies and the behavior of C in topologically nontrivial backgrounds would be wrong. 
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The mechanism by which lattice regularizations avoid contradiction with the known continuum 
properties can be seen already for the free field case and goes under the name of "doubling" [7]. 
Indeed, in four dimensions, relations (4.9), (4.20) and (4.21) require C to be of the form 

C{xaj,yPk) = Y,<^,i<^,PMU„Ul;fi)]{xj,yk) (7.2) 

with 

{[f{U„Ul;fi)]{xj,yk)r=[f{U;,U^;fi)]{xj,yk) (7.3) 

For U^ = 1 the simplest choice for / is 

[/(i> 1; M)](2;i, yk) = ^^jk[Sy,x+ii - Sx,y+ii] (7.4) 

Upon Fourier transformation of (7.4), C(p) becomes proportional to X^„ "'^/(p; m) with f{p;fJ.) = 
sinpi^ = pi_i. C(p) has 16 zeros rather than only one at p = as we would have liked. Linearizing 
around the zeros we get 16 copies of the continuum a ■ p but with some of the signs of the (t^'s 
switched. 8 of the (t^ sets thus obtained transform as the original set and the other 8 transform 
as the set aj^. We see that we have in reality not one Weyl fermion of a given handedness but 8 
right handed ones and 8 left handed ones. It is impossible to couple the link variables in a local 
and gauge invariant way to only some of these particles since they are represented by the Fourier 
modes of the same local field. Therefore, upon gauging one ends up with a vector theory and there 
is no contradiction with the know continuum properties. The structure of the regular lattice, as 
reflected in the toroidal topology of the momentum space (in the infinite volume limit) ensures 
that no other choice for the functions /(p; /x), even giving up some of the symmetries we required, 
will eliminate the doubling (probably a better term would be mirroring) of the Weyl fermions. In 
summary, C is not an acceptable discretization of the corresponding continuum operator because 
it describes too many fermion species avoiding the reproduction of anomalies and because it cannot 
have an analytical index. 

The above considerations mean that one cannot let the matrices B stay constants like in 
the formal treatment. If we did, the derivations presented there hold here too, and the overlap is 
simply related to the determinant of C, but the latter does not describe the theory we want. Let 
us first see how to rescue the overlap in the free case, but without changing C. We need to break 
the tight relationship between the overlap and det C for momenta away from the origin. If 6+ is a 
negative constant and C ~ 0, the vacuum of H"*" contains filled single particle states distinguished 
by having a vanishing bottom half; similarly if B~ is a positive constant and C ~ 0, the vacuum of 
H~ contains filled single particle states distinguished by having a vanishing top half. Thus every 
single particle state filled in one of the vacua is orthogonal to every single particle state filled in 
the other vacuum and the zero of C propagates to a zero in the overlap. We want this to happen 
only for the zeros at p = 0. Therefore, we make the matrix B space dependent, introducing a 
derivative term, such that in Fourier space B~(p) be positive everywhere, but B"'"(p) be negative 
only near p = and positive near all the other zeros of C: 

B'^(a;, 2/) = - ^[2^^,^^ - 6y^^+j:^ - S^^y+j:^] ^m < m < 1 (7.5) 
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This shows that the well known Wilson mass term eliminates unwanted fermionic excitations for 
the overlap also. Note that we have not put in any dependence on the spin indices; this helps 
preserving the symmetry properties we discussed above. We have also ignored the group indices 
since they are immaterial in the free case. 

We still need to show that upon gauging indeed the continuum behavior in perturbation theory 
and in nontrivial topological backgrounds is reproduced. This is the main goal of the entire paper. 
Here we just wish to explain that the reasons given before for C being an incorrect discretization 
no longer apply. In order to preserve gauge invariance when the gauge fields are turned on the 
derivatives appearing in B^ and C must be covariantized; therefore B^ and C no longer commute. 
It becomes now impossible to write down an explicit formula for the overlap in terms of smooth 
functions of the matrices B^ and C: the derivation given for the continuum, formal problem, 
breaks down. The overlap is still implicitly defined, but it is no longer guaranteed that one can 
make it smooth in the link variables expanded around C/^ = 1 and gauge invariant at the same 
time. It is also clear that the possible gauge breaking is restricted to the phase. With the Wigner 
Brillouin choice the overlap will still obey the set of invariances discussed in section 4 and this 
provides useful restrictions as mentioned there. Thus anomalies can be reproduced. 

Also, the possibility of correctly reproducing the effects of nontrivial topological sectors is 
opened up. For this to work one must see that unlike the free Hamiltonians, the ones with gauge 
fields turned on do not have impenetrable (gauge field independent) gaps in their spectrum. Indeed, 
on a lattice the space of the link variables is connected (the topological difference between the 
continuum A^ representing an instanton say and A^ = is lost upon passage to the link variables 
on the embedded lattice) so it must be possible for an eigenvalue to go through zero as the gauge 
field is changed. This has to be true of at least one of the Hamiltonians (we shall see later that 
indeed it is true only of H"*" but not of H~ and this is just enough). Unequal filling levels of the 
two vacua are able to ensure the robust vanishing of the chiral determinant reflecting the known 
related continuum property. 

In summary the two Hamiltonians are given by (7.1) with 

C{xaj,y(3k) = ^ a,(a,/3)^[^,,.+^C/;^» - S,,y+^{Ul{y)r] (7.6) 

B^{xai,yPJ) = -Sap Y.i'^^^J^i " ^y,^+f^K'('') " ^^,y+f^iUl{y)y'] Tm = BTm (7.7) 

and define the two many body ground states with phases fixed by the Wigner Brillouin choice 
relative to the appropriate reference states (depending on topology). The two ground states can 
then be used to compute the overlap or fermionic expectation values as appropriate. Regularized 
versions of the symmetries discussed in section 4 hold. The two Hamiltonians are local and typically 
have a gap. Each one of the ground states should be local in the sense that sampling it with various 
observables should generate local functionals of the gauge fields. Nonlocal terms in the action can 
appear once the overlap of the two different ground states is constructed. 

The action is known to include nonlocal terms of a special kind if the theory is anomalous. 
It is therefore important to see in more detail how anomalies appear in the overlap. We are 
considering here perturbative anomalies and therefore the background field is assumed to have zero 
topological charge.* Let |i2± >u be many body states obtained by filling all the positive energy 



To be sure, the anomalies have to appear also when the background is topologically non-trivial. 
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states of H^[U) respectively, with arbitrary overall phases. The Wigner-Brillouin phase choice is 
implemented by the replacement 

The chiral determinant for a single righthanded fermion is the overlap of the two ± states on the 
righthand side of (10.15). Eliminating a gauge field independent constant, the effective action 
induced by a single righthanded fermion becomes 

r^((7) _ i< -^ ~ 1-^— >u u< R — \R+ >u u< R + 1-^+ >i /„ Q\ 

~ \i< R- \R- >u\ i< R- \R+ >i \u< R+\R+>i\ ^ ' ' 

Under a gauge transformation, 

U^^{x) = g\x)U,{x)g{x + fl), (7.10) 

the Hamiltonians transform as 

Il^{xai,yPJ;U') = G\xai-x'a'i')^^{x'a'i',y'(5'j'-U)G{y'(5'j',y(5j) (7.11) 

where 

G{xai- yPJ) = g''{x)S^ySap. (7.12) 

Therefore the single particle eigenfunctions at U and U^ are related by 

^(C/s) = G^{U)e"^^^^3) (7.13) 

and the induced relations on the many body states are 

\R±>u.= g^\R±>ue"^^'-^'''> (7.14) 

Q is defined by G and is the same for both ± states. When (7.14) is used in (7.9) the unknown 
phases ^±{U,g) drop out. In the continuum the anomalies completely determine the dependence 
of the action on coordinates along gauge orbits. This dependence is given by the gauged Wess- 
Zumino Lagrangian [13]. If our regularization is correct it must also provide a lattice version of this 
continuum Lagrangian. The lattice Wess-Zumino action is a functional of the link fields U^i^x) and 
the group valued variables g{x). (When g is infinitesimal the dependence on U^ and g is restricted 
by Zumino's consistency conditions.) Define 

Then the lattice form of the exponent of the gauged Wess-Zumino model is: 

v-iW)-T-iu) ^ W*_{U,g)W+{U,g) 

\W*_{U,g)W+{U,g)\ ^ ' > 
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To obtain the anomaly one expands the logarithm of the above pure phase around ^r = 1 to linear 
order. There are two additive contributions, one from the negative side and the other from the 
positive side. As we shall see later, in the continuum limit only the contribution from the positive 
side survives. 

The locality of the H comes at some price: If one sets up perturbation theory from a path 
integral representation of the overlap one has to deal with a continuous time parameter and the 
action has a theta function discontinuity. This induces some technical complications that are 
likely soluble*. In our previous work we avoided the ambiguities related to the theta function 
discontinuity by discretizing the fifth coordinate also (we use language appropriate to 4-D, but 
everything holds in 2-D also). In this framework one deals with two transfer matrices roughly 
equivalent to the exponents of our Hamiltonians. The structure of the transfer matrices is quite 
similar to the structure of the Hamiltonians, but strict locality is lost. To be sure, the expressions 
are still local in the sense that the couplings between far away sites are suppressed exponentially. 

For completeness we recall here the relevant formulae for the transfer matrices obtained in [4]: 
The many body problems are still written as in (5.2) but the new H s are given implicitly, in 
terms of their exponents. There is no loss in doing this because the single particle problems are 
solved by diagonalizing the exponents now, the single difference being that the boundary between 
filled and empty states in the vacuum is at eigenvalue equal to one. Using our previous expressions 
for C and B^ we now have: 

By (7.7) is is clear that the matrices 1 + B are positive definite for any < m < 1. Therefore 
the inverses can be expanded around unity and one sees that the expressions, while not local in 
the strict sense are satisfactorily local, as mentioned above.** Numerically, this promises fast 
convergence of the conjugate gradient algorithm used to work out the inverses, thus the penalty for 
not having strict locality isn't that harsh as one might fear at first sight. Still, it is easier to work 
with (7.1). 



8. Topological charge on the lattice. 

The topological charge to be associated with a set of link variables is defined so that the 
fermionic reaction is the same as in the continuum: It is given by the difference in the number of 
negative energy levels of H"*" less the number of negative energy levels of H~. Clearly this is a 
gauge invariant definition. A direct consequence of the connectedness of the space of f/^'s is the 
existence of backgrounds for which at least one of the Hamiltonians has a zero eigenvalue. In the 
context of "lattice topology" such gauge field configurations are usually referred to as "singular". 

If the number of sites is V the dimension of H^ is 2 2 T^N X 2 2 T^N and, for U^ = In there are 
exactly 2 2 T^N/2 negative (positive) eigenvalues for each sign. As a matter of fact, H~ has 2 2 T^N/2 
negative (positive) eigenvalues for any gauge background. To prove this, we note that it is true for 



* B. Blok, private communication [14] 
** In our previous work [4-6] we used the notation B for the quantities 1 + B in (7.8). 
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the free case and next show that H can never have a zero eigenvalue: Let us assume that such an 
eigenvalue existed and that the corresponding eigenvector is 



B-u+Cv =0 
C^ii-B-i;=0 



(8.1) 



H(m) = ( ct -B .. ) (8.2) 



Multiply the first equation by u^ from the left and the second by v^ from the left. Complex conjugate 
the second equation, use the fact that B~ is hermitian and subtract the two equations. You get 
v^'B~v + u^'B~u = which contradicts the rather obvious positivity of B~. Similarly one can 
show that for any gauge field U^ all eigenvalues E of H~ satisfy u^[B~ — E)u+ ■!;^(B~ + E)v = 0, 
implying \E\ > m, so H~ always has a gap around zero. Clearly the proof breaks down for H"*" 
because B"*" isn't necessarily positive. 

Consider now B~ say but let m = /x be a varying parameter with the gauge background fixed. 
Denote the corresponding H~ by H(/x). As long as /x is positive no level has crossed zero; when 
/x becomes negative the gap can close. When /x = — m B~ becomes B"*" and H(/x) becomes 
H"*" . The imbalance of negative energy eigenvalues of H"*" is equal to the number of odd eigenvalue 
crossings occuring for H(/x) from positive to negative less the number of odd eigenvalue crossings 
from negative to positive when /x is monotonically varied from m to — m. A crossing is "odd" if the 
first nonvanishing derivative with respect to /x of the respective eigenvalue at the point where the 
eigenvalue vanishes is of odd order. 

Generically, crossings are odd with the first derivative nonzero. So, in practice, in order to 
determine what the topological charge of a given lattice gauge configuration all we need to do is to 
follow the eigenvalue flow of 

B + fi C 

through zero. By looking at the lowest eigenstate of H^(/x) and checking the expectation value of 
H/x) in the respective eigenvector a relatively efficient numerical method can be devised to measure 
the topological charge of a given gauge field configuration. 

Linear crossings of zero for eigenvalues of H(/x) occur at values /x where 

H(/x)V' = 0; VSV' 7^ (8.3) 

Rather than look at the hermitian eigenvalue problem (8.3) one can make the /x dependence explicit 
by considering the equivalent eigenvalue problem 

75H(/x)V' = => f ^^ ^')v = DV' = -/xV' V'SV'7^0 (8.4) 

We shall assume that all crossings are generic which means that the first derivative of the 
crossing eigenvalue is not zero at the crossing point. Then the last equation means that the 
topological charge can be obtained by looking at the usual lattice Dirac operator with a Wilson term 
(its coefficient, usually denoted by r, is unity - a preferred value for other reasons) and adding or 
subtracting unity for each reaZ eigenvalue in the interval [—m,m) depending on the sign oi ip^jsip. 
The lattice Dirac operator D is hermitian under the indefinite inner product (V'j^) = V'^75^- 
Therefore, D can have non-real eigenvalues A only if (^i/j\,i/j\) = 0. Thus, in the generic case, the 
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original definition in terms of the imbalance in the number of positive versus negative eigenvalues 
of H"*", is equivalent to: 

ntop = lim y^ —-r-^ — ; Btpi = Xitpi (8.5) 



{i; |3iA.|<m} 



IV'/tsV'^ 



D can be thought of as a lattice version of an approximately massless and antihermitian 
continuum Dirac operator. The above equations can then be viewed as a regularized version of the 
following continuum relation:* 

rztop = limV }l^^f ■ j-{d + iA)f, = \J, (8.6) 

This is another indication that the lattice definition of the topological charge obtained from the 
overlap is reasonable.** 

Clearly the definition contains some arbitrariness associated with the choice of m. In the 
continuum, this arbitrariness should have no effect if one measures the topological susceptibility for 
example. Most crossings reflecting real continuum instanton-like fields should take place for small 
values of the parameter /x. It may appear that we have an additional example of an ambiguity in 
our definition because we could use the transfer matrices instead of the Hamiltonians. However, a 
closer look reveals that the two definitions are identical if the parameter m here is identified with 
the parameter m there. To see this we need to go into some detail: 

First we show that, just like H~, T~ has an impenetrable gap around unity. Let us represent 

the eigenvector of T~ corresponding to the eigenvalue A by I ] . The eigenvalue equations 

-u Cv =\u 



1 + B- 1 + B- ^3_^) 

C^ u - C^ Cv - (1 + B-)v = -\v 

1 + B- 1 + B- ^ ^ 

are equivalent to: 

nt , 1 + B- 

(8.8) 



:tu: 


= - 


1 + B- 
"+ A 


— V 


Cv: 


=u 


-A(l + B- 


)u 



* In (8.3) (8.4) and (8.5) the notation ip^jsip implies also a sum over the sites x, while in (8.6), 
following continuum conventions, the integral over x is made explicit. 

** Another fermionic definition of the topological charge has been given in [15]. These authors 
elected to focus on the continuum formula ritop = m J trlj^ ./-gi^^-)!^ ] replacing it by ritop = 

i^p'm jY^i m+\-M, where D/^ = \ifi, /^D = \ifi and / fifi = 1. For the full definition of the 
parameters Kp and Mc we refer to the source. For smooth backgrounds one can take Kp = 1 and 
Mc = —miriii^^Xi). The parameter m is to be extrapolated to zero. For such smooth backgrounds 
we expect good agreement with our definition. 
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Multiply the first equation by v^ from the left, the second equation by u^ from the left and complex 
conjugate it. Subtracting the so obtained equations one obtains: 



,t 



1 + B- 
A 



v + u^ [A(l + B-)- 1] 1^ = 



(8.9) 



The above is impossible if A is such that B + 1 > max[\, ^). Since B > m there is a gap for 

1 



1 + m 



< A < 1 + m 



(8.10) 



Thus, just like in the Hamiltonian case, the topological charge in the transfer matrix formu- 
lation will again be given by the number of (assumed generic) crossings of unity of eigenvalues of 
r(/x) as we vary /x from — m to m. ^(/x) is defined as T~ with the parameter m replaced by /x, 
so /x interpolates between T~ and T"*". We shall show now that generic crossings occur for H(/x) 
and ^(/x) at the same time and with the same sign, and thus obtain the equivalence of the two 
definitions for the topological charge for generic gauge configurations. 



The unit eigenvalue condition for ^(/x) where the eigenvector is taken as i/j 



IS 



C^ix-B-i;=0 
Ci; + B"ix=0 



(8.11) 



This equation is identical to (8.1). To see the relation between the signs we need to compute the /x 
derivative of the crossing eigenvalue at the crossing point. Using the eigenvalue equation it is easy 
to see that 
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(8.12) 



which is the desired relation. This concludes the derivation of the equivalence of the two definitions 
of the topological charge. The continuum interpretation of the hamiltonian definition thus holds 
for the transfer matrix definition also. 

For H(/x) it was useful to consider H^(/x) because for it the states which will potentially cross 
are located near the lower end of the spectrum. When working with the transfer matrix version 
of the overlap the role of H^(/x) is played by ^(/x) + r~^(/x). It is easy to find the inverse of the 
transfer matrices in (7.8): 



(r±)-i=exp(-H±) 



Ci+^Ct + M 

l + Bi*-^ 
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(8.13) 



We see that the inverses are just as local as the transfer matrices were; this comes as no surprise 
since the inverses generate propagation in the opposite direction in the auxiliary "time" space. 
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9. Construction of anomaly free theories 

The ground work in the previous sections makes the procedure for constructing vector and 
chiral theories quite straightforward: Sections 2 and 4 have estabHshed the formaHsm for treating 
a single right (left) handed fermion coupled to a gauge field on the lattice. The method to compute 
fermionic expectation values was described in section 5. The lattice regularized theory was described 
in section 7. This forms the basis for constructing vector theories and chiral theories. In this 
section we construct such theories and show that several fundamental properties of these theories 
are reproduced. 

We start by constructing a single flavor massless vector theory. We decompose the Dirac 
fermions in the single gauge multiplet into their left handed and right handed components. To each 
chiral part [R or L) we associate a set of creation and annihilation operators [a;^„j] ; [c^l^aj] i ^x^ai'i 
a^ai- Both sets separately obey the canonical commutation relations given in (5.3) and the two sets 
anticommute with each other. From each set of creation and annihilation operators we construct 
a many body Hamiltonian as given in (5.2). The single particle Hamiltonians appearing in both 
of these many body Hamiltonians are the same and are given by (7.1), (7.6) and (7.7). The single 
particle Hamiltonians are isomorphic because both the right and left handed fermions are in the 
same representation of the gauge group, carry the same charge and are coupled to the same gauge 
field. The two ground states that constitute the overlap are denoted* by |F± >u. |F± >u are 
direct products of |i2± >u and |X± >u where |i2± >u (|-Zv± >u) are the specific many body states 
associated with the Hamiltonian of the right (left) handed fermion as described in the beginning 
of section 5. The overlap is a product of the two overlaps, one for the right handed and the other 
for the left handed fermion and are given by (5.4). 

To compute expectation values of fermions we first decompose them into left and right compo- 
nents. In the massless case the left and right contribution factorize and the individual expectation 
values are then computed by following the procedure described in section 5. The generating func- 
tional is given by 

Z^(7?^, 7?^, 7?^, 7?^; U) = u<V-\ exp [f/^a^ + 7?^[a^]t + 7?^a^ + vV]^] \V+ >u (9.1) 

In the above equation, we have suppressed the sum over xai explicitly shown in (5.10). In the zero 
topological sector the result is of the form 

Z^'ifj^, r,^, f,^, r,^;U) = u<V-\V+>u exp [-7?^G^(C/)7?^ - 7?^G^(C/)7?^] (9.2) 

where G^{U) = — [G^(C/)] and the explicit expression for G^{U) is given in (5.22). It should be 
kept in mind here that the sources are for Dirac fermions and that only half the modes of G^[U) 
propagate and these correspond to the physical right handed chiral fermion. 

The massive vector theory is thought of as a perturbation of the massless case. The mass term 
that couples the right handed fermion to the left handed fermion is treated in a manner similar to 



* The gauge fields are denoted by the link variables U to show that the definition is on the 
lattice. 
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the fermionic sources. The generating functional for a massive theory is written as 



Z^{fi^, r]^, 7?^, 7?^; U, mf) = u<V-\ exp mf[a^]^a^ 



ta"^ + 7?^[a^]t 



exp 
exp 



7] a + 7/ [a J ' 



mjla^^^a^ 



|F+> 



(9.3) 



u 



where mj is the fermion mass (here assumed to be positive for simplicity). The mass term can be 
rewritten in terms of fermionic source terms: 



Z^(7?«,7?«,7?^7?^; U,mj) = / [de][de][di%di%< F - l exp V^(^"[a«]t + i^a 
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(9.4) 



In the zero topological sector (9.2) can be used in (9.4) and the integration in (9.4) can be 
carried out. This results in the following expression for the generating functional: 



Z^(7?^,7?^,7?^,7?^;C/,mj) = u<V - |F+ >u det if (C/,mj) exp -fj" G" {U,mj)r] 



where 
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The matrix K{U,mj) in (9.5) is given by 

K{U,mf) = 
where 



mjG^U) 



-mjG^U) 



From (9.5), (9.7) and (9.8), the expression for the vector determinant is 



(9.7) 



(9.8) 



j<V - |F+ >u det K{U, mj) = u< R - \R+ >u u< L - \L+ >i 



det 



L^T-rMt 



l + m}G'^{U)[G'^{U)] 



(9.9) 



Since u< R — \R+ >u= u< L + \L— >u (c.f. Lemma 4.1) it follows that the above formula for 
the vector determinant is real and positive for all gauge fields. This is true at the level of a single 
flavor and holds whether the fermion is massless or massive. Since the determinant is real it is also 
gauge invariant because the only gauge breaking in the right and left overlaps were in their phases. 
(9.9) should be compared with the formal expression for the vector determinant 



det D = det C det C^ det 



l + mlC-^[C 



tl-i 



(9.10) 
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The matrix G^ {U^mj) in (9.5) is given by 



'^'•<--)=(-J(.,rT) 



\-mi[G^{U)\ . r. nt "T'/ — r. It- / 

\ -^ l + m'jG^iU)[G^iU)\ l + m'j[G^iU)\ G^iU) / 



(9.11) 



The first term is the propagator for massless fermions and the second term includes all the effects 
due to the mass. In a manner similar to the relation between (9.9) and (9.10), the Dirac propagator 
(9.11) also has a formal analogue expression. The diagonal terms in (9.11) lead to a version of the 
well known [16] spectral criterion for ^^ to have non-zero expectation value when rrif -^ 0: 
Let the eigenvalue density of G^[U)[G^[U)] , after averaging over gauge fields (including the 
fermion determinant), be /o(A) where the eigenvalues A are non-negative. ^^ maintains a nonzero 
expectation value as rrif ^ if for large A we have /o(A) ~ tjj^- 

Turning now to an arbitrary number of flavors, Nf, amounts to a trivial extension. We can 
think about the physical Dirac fermions as basic objects whose propagator is G^ {U,mf).* The 
extension of relation (9.8) (equivalently (5.13)) to the Dirac propagator is 

75G^(C/,mj)75]^ = G^(C/,mj) (9.12) 

The above equation is satisfied by (9.11) because the diagonal terms are hermitian and the off 
diagonal terms are related by one being the negative hermitian conjugate of the other. Having this 
equation hold in a regularized framework independently of whether the theory is massless or not 
is useful when more flavors are considered because positivity of the regularized determinant (c.f. 
(9.9)) and relation (9.12) are the only ingredients needed to prove the well known mass inequalities 
between mesons and baryons. The original proofs in the path integral formalism go through directly 
[17]. In the lattice case, the improvement relative to Weingarten's work is that a restriction on the 
number of flavors is no longer needed.** 

The massless multiflavor case has at the formal path integral level a U{Nf) X U{Nf) symmetry, 
the two factors acting on the left and right components independently. It is important to see 
how these symmetries are realized after regularization by our method. This can be easily done 
because the symmetry is there for each gauge field, even before carrying out the full functional 



* One may easily replace the common parameter rrif by a more general mass matrix connecting 
each physical righthanded gauge multiplet to some identical multiplet of lefthanded nature. This 
matrix can be brought into diagonal form with the help of some canonical transformations on the 
creation/annihilation operators. 

** Traditional lattice regularizations of vector theories involve a parameter k that has to be fine 
tuned to a specific value Kc (dependent on the lattice gauge coupling) in order to have finite or 
zero mass quarks in the continuum: While the numerical values of Kc are not known exactly, we 
know that when the continuum limit is approached and the lattice gauge coupling goes to zero, Kc 
is such that the positivity of the fermionic determinant per individual flavor is lost and one needs 
an even number of flavors to preserve the needed positivity of the gauge integration measure. 
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integral. The symmetry can be seen in the formulae for the fermion generating functional (5.9) when 
properly generalized to many flavors. For each one of the symmetry transformations an associated 
canonical transformation on the creation/annihilation operators a, a^ exists and is independent of 
the background gauge fields. If the total vacua |F± >u are both annihilated by the generator 
of the above canonical transformation the generating functional will exhibit the symmetry by its 
functional dependence on the sources fj,rj. For the SU{Nf) X SU{Nf) part of the symmetry group 
it is easy to ascertain that the symmetry indeed holds. The unimodularity of the matrices in flavor 
space plays a role in deriving the invariance of the two vacua. Similarly, one can easily convince 
oneself that the vectorial U{l)v is also preserved, although the two vacua now rotate by a phase. 
Indeed, the phase coming from \V— >u exactly cancels the phase coming from \V-\- >u and overall 
invariance holds. This is true even if there are individual imbalances in the level of filling due to 
nontrivial gauge backgrounds. The effects of the imbalances in the righthanded and lefthanded 
sectors cancel each other. For an axial U{1)a transformation the two imbalances add to each other 
and therefore U{1)a is not a symmetry if the background has topological charge. This is just what 
we would have liked to see happen. Let us take a background carrying topological charge equal to 
unity. The axial transformation formally known as ^ ^ g«"75i^ translates into a^'^ -^ e~^"a^'^ 
and a^'^ -^ e'"a^'^ for / = l,2,---iVj. The vacua transform as \V+ >u^ e'^'"^f\V+ >u and 
\V— >u^ \V— >u- For a fermion expectation value to be nonzero it has to transform in the 
opposite way. That is what a 't Hooft vertex does. 

If Nf = 1 the 't Hooft vertex is, in terms of physical fields, given by ^'^^l, which together 
with its hermitian conjugate becomes the mass operator. If there are Nf flavors, the 't Hooft 
operator is made up of 2Nf fermi fields. When the 't Hooft operator aquires an expectation valus 
this indicates a breaking of the axial U{1)a, down to a discrete subgroup. The breaking occurs 
even when the volume is finite and is clearly distinct from a spontaneous breaking. The global 
symmetry Ur{^N j) X Ul{^N j) formally present in the massless case gets explicitly broken down to 
SUr{Nj) X SUl{Nj) X Uv{l) X Z2. 

When a mass term is turned on terms with the quantum numbers needed to give a nonzero 
matrix element between the two vacua for any gauge background will be generated by expanding in 
the mass. Thus, the vectorial determinant will not vanish for any topological charge, as expected. 

We have shown that our regularization preserves the topological effects one is familiar with from 
the continuum as far as quantum numbers go and we shall see more of this below. However, we ought 
to remember that the continuum reasoning is credible because one can imagine setups where one can 
carry out semiclassical computations with confidence and indeed show that quantities that could be 
non-zero indeed are. (For example, we may enclose the system in a small Euclidean volume, so that 
the effective couplings never get strong.) The magnitude of some of the above topological effects 
can then be computed using instanton techniques. If our regularization is indeed right we have to 
convince ourselves that one could carry out the instanton computations directly on the lattice and 
that the results would be the ones we know from the continuum. The main point to check is whether 
the overlap indeed behaves in a way similar to the continuum determinant when the background is 
a "dilute gas of instantons and anti-instantons" [18]. The dilute gas approximation was important 
in at least two respects: It showed explicitly how clustering is restored after integration over all 
topological sectors and by this reassured us that indeed a unitary local field theory can incorporate 
the somewhat unfamiliar situation of having to carry out a functional integral over a disconnected 
set of fields of the same type. The dilute gas approximation also indicated how a theory with no 
massless particles, like QCD for example, while having no sensitivity to the boundary conditions still 
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manages to exhibit effects refiecting topological properties of the gauge field space. The dilute gas 
approximation achieved the above by reintroducing extensivity in the thermodynamical sense via a 
factorization of the chiral determinant into approximately local factors when the gauge background 
consists of widely separated instantons and anti-instantons. Similar issues arise when one considers 
below the insertion of 't Hooft vertices and their exponentiation. In the following sections we shall 
try to verify that indeed the overlap behaves as a determinant in the sense that it factorizes in 
dilute backgrounds. We shall also try to see how effects depending on the explicit form of the 
continuum fermionic zero modes get reproduced by the overlap. 

The definition of the phases of the states given in section 5.2 for all topological sectors also 
applies to the vector theory with N j flavors. J^ used to define the phase is the 't Hooft operator that 
breaks the axial Ua{X) symmetry. In the definition T can be replaced by e^^ T with 6 being some 
real parameter. This just replaces the reality condition of the Wigner-Brillouin in non-zero sectors 
with a condition where the phase is fixed at nO in the sector with charge n. The only difference 
would be that different topological sectors will now be weighted by e'"^ which is equivalent, in 
four dimensions, to the addition of a -^^ Jtr[FF) term to the -^ J tr[FF) term of the gauge 
action. This is the usual 9 term that one would like to introduce in a vector theory. The phase 

factor in front of J^ can be absorbed by redefining the fermion operators by a^'^ -^ e ^"f a^'^ 

and a^'^ -^ e ^"f a^'^ with X^jCC/ = Nf. The many body Hamiltonians are invariant under this 
redefinition and only the mass terms get affected. The ability to move the ^-parameter dependence 
from the J FF term to the mass matrix is familiar from the continuum and we see that it holds 
exactly in our regularization. In particular, by appropriately choosing the constants cc/ we see that 
the vacuum energy will be 9 independent if one of the Dirac fermions is massless. 

The construction of a chiral theory follows in the same manner. Each multiplet of a given 
handedness is represented by one of the building blocks we have constructed. Since the right and 
left handed fermions are not paired in a chiral theory the resulting overlap need not be real whether 
the theory under consideration is anomaly free or not. As will be verified later on, each block carries 
the correct amount of "consistent" anomaly. The difference between an anomaly free theory and 
an anomalous one is that the leading operators in the imaginary gauge breaking terms cancel in 
the former (i.e. the coefficient of the induced gauge Wess-Zumino action vanishes in the anomaly 
free case). However, even in the anomaly free case exact gauge invariance does not get restored for 
finite lattice spacings: Higher order terms that break the gauge invariance could still be present on 
the lattice, but their effects are expected to disappear in the continuum limit. This point will be 
discussed again in a later section. 

The standard model can be constructed by following the same procedure: The new ingredient 
are the Yukawa interactions and they can be added in a manner that is similar to the addition of 
the mass term in the vector case. There is one major difference beyond the obvious space time 
dependence and gauge group representation: The mass term in the vectorial case would eliminate 
all the zeros of the fermionic determinant in backgrounds carrying non-trivial topological charge. 
However, in the standard model the Yukawa term couples fermions transforming as doublets under 
SU(2) to fermions transforming as singlets under SU(2), so only one of them sees the topological 
charge of the SU(2) gauge field. Therefore, no term in the expansion in the Yukawa interaction 
can undo the vanishing of the overlap. One still needs to insert a 't Hooft vertex to get a nonzero 
matrix element. This feature is important in understanding baryon decay in the standard model. 
There have been several recent papers proving the resilience of the topology induced zeros when 
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there are Yukawa interactions present [19]. 

If we consider SU{N) gauge theory interacting with a single adjoint multiplet of lefthanded 
Weyl fermions (gluinos) the continuum limit should come out to be supersymmteric. Unlike in 
previous approaches no fine tuning is needed because the fermions are strictly massless (see [20]). 
However, we have not yet detected any remainders of the continuum supersymmetry on the lat- 
tice. Still, we would hope that this observation is an indication that progress on the problem of 
constructing lattice models that have some remnants of supersymmtery could be made using our 
methods as an inspiration. 



10. Abelian Gauge models in Two Dimensions 

Massless fermions in two dimensions differ in several related respects from massless fermions 
in four dimensions: Anomalies can be cancelled only between left and right movers. Charge conju- 
gation does not turn a left mover into a right mover. Anomalous gauge theories are renormalizable. 
Anomaly free theories can be regularized to all orders in perturbation theory gauge invariantly by 
adding a finite number of Pauli-Villars Dirac fields of wrong statistics.* Abelian models are soluble 
since the determinants are exactly known, and moreover, turn out to be quadratic in the gauge 
fields rendering all the models of massless fermions interacting with abelian gauge fields gaussian 
in the trivial topological sector. Like in four dimensions, some information on the slightly more 
complicated higher sectors can be obtained from the zero topological sector by assuming clustering. 
Even truly chiral abelian gauge models (with nontrivial anomaly cancellations - not of the vector 
type) have no parity violating terms in the chiral determinant, i.e., in Euclidean space the action 
has no imaginary part. 

In spite of the above we believe that tests of the overlap regularization in two dimensional 
chiral abelian models are of significance. The objectives of the tests fall into two classes: One is to 
establish that the lattice overlap correctly reproduces the chiral determinant (both its imaginary 
part and its real part) and the salient features of the fermion two-point function in backgrounds 
that carry zero total topological charge. This will be the subject of subsections 10.2 and 10.4. The 
second is to show that the lattice overlap correctly reproduces the continuum features of fermions 
in topologically non-trivial backgrounds. This has to be done at the "classical" level, i.e. showing 
that the correct effects of fermionic zero modes with the right quantum numbers are reproduced 
(in subsection 10.3), and also at the "quantum" level. By "quantum" we mean the needed features 
that enabled 't Hooft to replace the fermionic effects caused by instantons by an effective local 
vertex. To very good accuracy, in four dimensional spontaneously broken gauge theories, if 't 
Hooft 's vertex is included instantons can be ignored, and the theory is perfectly "normal", obeying 
clustering. The replacement works not only because one has the right number of zero modes, but 
also because of the form of their asymptotic decay at large Euclidean distances from the instanton 
core. Thus, it is important to show, at least in some simple cases, that some general properties of 
the shape of the continuum zero modes are preserved by our lattice regularization. This will be 
done in subsections 10.4 - 10.6. In subsection 10.7 we shall show that non-trivial global topology, 
and the associated zero modes, do not spoil the perturbative anomaly equations. 



We thank S. Yankielowicz for pointing this out to us [21]. 
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The simple tests that we have carried out for specific backgrounds strongly indicate that the 
complete dynamical regularized quantum model will have the correct continuum limit. We know 
of no other regularization method that has gone as far. We believe that any credible method for 
latticizing chiral fermions with gauge interactions is required to have an explicit mechanism showing 
how continuum instanton effects are reproduced. The paragraph at the end of section 6. shows by 
example how even seemingly perfect perturbative behavior is still an insufficient test. At the level 
of quantum numbers this requirement has been emphasized in [22]. 

The tests we shall present in detail below will all be for the Hamiltonian, local version of the 
lattice overlap. Similar tests have been carried out in some cases in the transfer matrix formalism 
and were equally successful. In particular, results of tests of the anomaly and eigenvalue flows dealt 
with in subsections 10.2 and 10.3 below have been previously reported in [5] and [6]. We shall always 
work on tori and use antiperiodic boundary conditions for the fermions. These conditions will be 
implemented by adding the appropriate U{1) gauge factors to the link variables. These phases will 
never be shown explicitly in our equations; whenever a link configuration is constructed, the fields 
displayed in the equations differ from the ones in our programs by the above phase factors. 

10.1 Review of Chiral Schwinger Model [23], [21] 

The partition function of the Chiral Schwinger Model is formally given by: 



dA^dipdipe 



Jd^4T^^'^-^-Er=i''^^<^^^+"i^^^>^^-ELi''^'^^^^^+"i'^^>^'^ 



(10.1) 



(Ti = 1 and (72 = i as defined in (2.6). We have included several single component gauge invariant 
fermions of both chiralities. We are tacitly assuming that we are working on a torus and that 
the infinite volume is to be taken whenever possible. The parameter g has dimensions of mass 
and so does the vector potential A^. The field F^j, = 5^Aj, — 5j,A^ = e^j,£^ has dimensions of 
mass squared and the topological invariant ^ J d^xE takes integral values. The fermion charges 
are integers also. The path integral over A^ includes a discrete sum over all the sectors. In the 
following we shall focus on expectation values of operators to which only the zero topological sector 
contributes. This simplifies the evaluation of the chiral determinant. 

Gauge fields of zero topological charge can be assumed to have a field strength that vanishes 
at infinity. The gauge potential can be decomposed into two real dimensionless scalar fields that 
also vanish at infinity and are defined by: 

A^ = d^X + e^.d,(j) (10.2) 

The integral over A^ is replaced by an integral over x a^nd (j). x should be thought of as an angular 
variables and if the theory were gauge invariant, the integral over x could be dropped for gauge 
invariant observables. There is no field dependence in the Jacobian associated with the A ^ X:4' 
change of variables. The following change of variables is done for the fermions: 

^f = e-'«?Xggf 0^/?' ^f = e»9f Xg-gf 0^/?' ^f = g-'^f Xg-gf 0^^' ^^ = g'^f Xg^f 0^^' (iq.S) 
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Under these changes of variables the action changes to: 



d^x 



. ^ r=l 1=1 



(10.4) 



To complete the change of variables one needs the Jacobian of the transformation on the fermions. 
Naively it would be unity, but Lorentz invariant regularizations will spoil this. 

The Jacobian comes from the ratios of the appropriate chiral determinants and their free coun- 
terparts. Since we are in zero topological sector we can evaluate the determinants in perturbation 
theory. Note that the right-handed determinants depend only on a* ■ A = Ai — iA2 while the 
left handed ones only on a ■ A = Ai + iA2. Any one loop fermion diagram with three external 
gauge fields or more is convergent and therefore (we assume that a regulator has been introduced 
and subsequently removed) unchanged by regularization. In particular it is gauge invariant. But 
one cannot make a gauge invariant quantity out of Ai — iA2 alone, without using also Ai + iA2. 
So all orders higher than the second vanish. We cannot demand the regulator to preserve gauge 
invariance and make the two loop diagrams finite for arbitrary charges. But we can find regulators 
that ensure the gauge invariance of the absolute value of the chiral determinants. Also, a regulator 
can be chosen to make the absolute value of the determinant parity even and the phase parity odd 
as described in section 4. It goes without saying that the regulator is chosen to preserve Lorentz 
invariance. These requirements when combined with dimensional considerations limit the form of 
the logarithm of the chiral determinant to: T^ | [A] = J d?x[ciE-^E + ic2E-^d ■ A\ where ci^2 are 
real constants. Their values are determined by the anomaly equations: 



?2 _^ T . r 2 
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(10.5) 



Note that the regulator induced dependences on both Ai ± iA2 in the chiral determinants: it is 
easy to see that a Pauli Villars regulator, for example, will do this. 
Therefore the Jacobian is given by: 



J = exp 



E. Q? + El <lf f ^2^^g2^ ^ -ErOl ElQl_ , ^2^02^ 
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(10.6) 



Our decision to choose a regulator that puts all the gauge breaking into the parity odd part of the 
chiral determinant (i.e. in its phase) fixes the parameter a of Jackiw and Rajaraman [23] to unity 
for unit charge. The anomaly free case where ^^ q^ = ^^ qf" = Q^ is a gauge invariant theory 
and has several additional U{1) symmetries, depending on the number of fermions. One of these 
extra global f/(l)'s is always anomalous. 

The most famous case is the original Schwinger model where R = L = 1 and q^ = q^ = 
1. Here the anomalous U{1) is axial charge. The symmetry is explicitly violated in nonzero 
topological backgrounds which give a non-zero expectation value to the axially charged, gauge 
invariant, Lorentz scalar operators V = 'ip^ip^ and V = tp^ip^. This can be seen by computing 
the expectation value < V{x)V(^0) > which gets contributions only from the zero topological 
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charge. Using the transformation to free fields one establishes that when cc^ ^ oo the correlator 
< V{x)V(^0) > does not vanish, implying via clustering that < V >^ 0. Since in the zero topological 
sector < V >= 0, clustering requires the inclusion of all topological sectors in the path integral. 
The calculation in more detail is: 



< V{x)V{0) > = < i;^\x)il;^\0)il;^\x)i;^\0) >< e2(0(-)-0(o)) 
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The integral in the last line decreases exponentially as cc^ ^ qq and the constant 7 is .577215.... 
From this one gets, assuming clustering, the well-known result: 



<V>\ = ^(^) 



(10.9) 



The simplest example that contains three fermion fields will have two globally conserved f/(l)'s 
and another global U{1) will be anomalous. The anomalous U{1) can be chosen to be the fermion 
number. To preserve local gauge invariance one needs a Pythagorean relation between the charges. 
The standard example is the 3-4-5 model: R = 2,L = 1, q^ = 3,q^ = 4:,q^ = 5,Q = 5. An 
analogue of the operator V of the Schwinger model is F = V-^'V-^V^: 
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(10.10) 



These operators are gauge invariant in spite of employing ordinary rather than covariant derivatives, 
because the extra terms vanish by antisymmetry. Under an Euclidean "boost" ip^ -^ e^^'^ip^, 



ijjR -^ e'"/2,/,-R 



Tp^, Tp^ 



-ia/2„l,L 



lP^,lP^ 



-ia/2„l,L 



ip^ anda-d ^ e'"a ■d,a*-d^ e-'"a* ■ d. Therefore 



F is a Lorentz scalar. V is obtained from V by interchanging every i/j with the corresponding i/j 
and every a-d with a* -d. V is gauge invariant and also chargeless under the global non-anomalous 
U{1) which acts by 
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e^'^^^, 



i'i 



R 



-itiq- 
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V^: 



R 



-lUq^ 



,j«gi 



Vr 
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(10.11) 



^L _^ ^L^ ^L _^ ^1 



When substituting the fermion fields in V and V with the free fields and the (j) exponentials 
the derivatives have to be taken as acting on the free fermion fields, to avoid vanishing by Fermi 
statistics. Denoting by V and by V' the expressions one obtains from V and V when one substitutes 
the ■^"s for the corresponding ■^'s one finds: 



< V{x)V{0) > = < V'{x)V'{0) > e 



4Q2g(! 



'-) 



(10.12) 
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Free field theory leads to: 



< V'{x)V'{0) > 



(27r)l«fl+l«2'l+l«^l(a;2)Q2 

kl-i 
/(|^|)= det Ki + jY] =111 (''■)]' 

0<l,'i<\q\—l ■^■^ 



(10.13) 



Exactly like in the Schwinger model one determines from clustering and the large x^ asymptotics 
of G that < V >^ 0. Since V carries fermion number 2 this means that the model will allow for 
reactions in which the total fermion number changes by two units. These reactions will take place 
in backgrounds that carry topological charge 1. 



' < ^ > ' - y (27r)kfl+k.^l+k^l ['^J ^^ ' 

For completeness, let us record the effect of finite temperatures on the expectations | <V> \. 
The derivations are essentially identical and one ends up having to replace the factors Yt" ^^ 

zero temperature by ^}^%- exp[jQ c?t(1 — e^Tvr'^"'* ^)~^] at finite temperature T. For very high 

temperatures the new factor simplifies to 2'KTe oQ . The finite temperature results are important 
in two respects: Simulations of the full dynamical system at finite temperatures would be easier, 
since the number of lattice points needed might be quite small. In addition, the behavior of the 
chiral condensate of the ordinary Schwinger model, < V + V >, with temperature emphasizes the 
fact that the axial symmetry is broken explicitly by the anomaly and not spontaneously (at finite 
temperature the system is effectively one-dimensional and any spontaneous breaking of a symmetry 
that occurs at zero temperature would have had to completely disappear). 

Both in the Schwinger model and in the 3-4-5 model and both at zero and non-zero temper- 
atures the U{1) violating processes can be understood as a consequence of the zero modes of the 
various chiral Dirac operators in a gauge background of nontrivial topological charge. Since the 
integral over the scalar field is Gaussian, for each observable there is one particular background 
gauge configuration that contains all the dependence on the observable. The expectation value 
of the observable is determined by one particular gauge field configuration and the fermionic de- 
terminant in it. If the chiral fermion regularization can be shown to preserve continuum features 
like the approximate zero mode structure in some fixed backgrounds made out of superpositions 
of instanton and anti-instanton like configurations, it is quite likely that the full model, including 
the fluctuations of the gauge field, has been regularized in a way that reproduces all the continuum 
features. We already pointed out that due to the quadratic nature of the induced gauge field action, 
once anomalies cancel, even in a chiral model, the chiral determinant is real, i.e. no parity violating 
processes between the gauge bosons are induced. Thus in the two dimensional abelian models, 
genuinely chiral models and massless vector models do not differ as much as they do in four dimen- 
sions. For this reason, any regularization that preserves the factorization of the chiral determinant 
into contributions from each fermion of each charge and chirality and has correct parity behavior 
would most likely reproduce the anomaly-free chiral models correctly if it reproduces the vector 
model correctly. In the following subsections we will show that the overlap regularization of the 
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abelian models reproduces all the necessary perturbative and non-perturbative features in fixed 
gauge backgrounds. 

10.2 Perturbative gauge backgrounds 

We intend to check that the lattice regularized overlap reproduces correctly the perturbative 
features. We shall pick a continuum smooth background on a torus given by a plane wave with a 
small number of nodes and a small amplitude. This background will be discretized using sets of 
lattices of decreasing coarseness. We wish to show that in the limit the overlap will be equal, as 
a complex number, to the continuum formulae reviewed in the previous subsection. Also, we shall 
show that when we vary the background by a gauge transformation the variation of the imaginary 
part is correctly reproduced, again in the continuum limit. (The real part is strictly gauge invariant 
and therefore unaffected by the variation.) It is sufficient to verify this for one chirality and one 
charge. We pick right handed fermions with unit charge. 

We embed an X X i lattice in the continuum torus whose physical size is I X I. The lattice 
spacing a is given hj aL = I. To compute the effective action in perturbative backgrounds, we 
choose 

A^{x) = Alzos^-j— (10.15) 

as our gauge field on the continuum torus. A^ has dimensions of mass and is held fixed. The k^ are 
fixed integers and determine the momentum of the gauge field. The electric field, E = diA2— d2Ai, 
associated with this gauge field is given by 

E{x) = -^{ka',-k2A',)sm^^!^^ (10.16) 

Referring to (10.6) we see that, ignoring field independent constants, the effective action for this 
particular gauge field is given by 



^'"^[^ ^^)] = - 8.{kl + kl) 

. (fciA;/ + fc2A^O(fciA^/-fc2A;o 
^"^"^[^ ^^)] = 8.{ki + ki) 

and the gauge variation of the effective action is given by 



(10.17) 



If we choose 



then (10.18) gives 



^ = -(,,4-feA;)sm^ (10.18) 



X{x) = x''sm^—- (10.19) 



^5r%4)_ 1 



(10.20) 



On the lattice, the link element corresponding to the gauge field (10.15) is 

r /-i ^ ] r.M°;sin^ p^ i r.M°; p^ 

U^[n) = exp i / A^[x + tafi)dt = exp i— ^ cos[p ■ n -\ ) ~ exp i—— cos(p • n -\ ) 

L Jo -' '- -^ ^ 2 J L X 2 

(10.21) 
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with lattice momentum 



2'Kk, 



Pi^ 



0<ki,« L. 



(10.22) 



n = x/a labels the discrete sites on the lattice. Note that the lA'i are dimensionless numbers. The 
gauge transformation (10.19) is represented on the lattice by 



Kn) 



Jx sin(p-ra) 



(10.23) 



We compute the overlap on several lattices with increasing L's and extrapolate to the contin- 
uum, L -^ 00, with the intention to verify (10.17). An exhaustive verification for all lA'i and A;^ 
would establish that the lattice regularized overlap reproduces all the exact continuum results for 
perturbative gauge fields. 




Figure 10.1 The real part of the induced action as a function of number of lattice points for a 
particular smooth configuration. 



As an example we choose lA^ = IA\ = 0.32 and ki = 1 and k2 = 0. We set m = 0.9 
in the Hamiltonians H . Recall that the parameter m is an additional ultraviolet cutoff in our 
regularization procedure. This parameter does not have to be fine tuned, except being kept large 
relative to physical scales. The continuum limit is independent of the precise value of m. The 
effective action is computed on a finite lattice by using (7.9). The finite dimensional Hamiltonians 
are diagonalized for the gauge field background given by (10.21) and the corresponding single 
particle eigenfunctions are obtained. From this the many body states |i2± >u are constructed. 
The real part is computed for various X's and plotted as a function of j- in Figure 10.1. Extraction 
of the continuum limit yields a value of —0.00406(2) which is in good agreement with the righthand 



side of the first line in (10.17), namely 



(0.32)^ 

Stt 



-0.0040744. The imaginary part is computed 



for various X's and plotted as a function of j- in Figure 10.2. Extraction of the continuum limit 
yields a value of 0.00409(2) which is in good agreement with the righthand side of the second line 
in (10.17), namely i^|^ = 0.0040744. 
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Figure 10.2 The imaginary part of the induced action as a function of number of lattice points 
for a particular smooth configuration. 



D -0.0005 




Figure 10.3 The contribution to the anomaly coming from the "inactive" side, A as a function 
of number of lattice points for a particular smooth configuration. 

The gauge variation of the effective action on an X X i lattice, following (7.16), is expressed as 



-K dx° 



x°=o 



A{U, L) = A- {U, L) + A+ {U, L) 



with 



'^'(^'^) = ^^^^°^ 



u 



< R + \g\R+ >i 



\u< R+ \Q\R+ >i 



(10.24) 



(10.25) 



A {U, L) are the contributions to the anomalies from |i2± > respectively. Q is as defined in (7.10)- 
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Figure 10.4 The contribution to the anomaly coming from the "active" side, A'^ as a function 
of number of lattice points for a particular smooth configuration. 

(7.14) with the g appearing in (7.10) given by (10.23). It has been observed before that H"*" is the 
"active" side. For instance, B"*" changes sign as a function of momentum in order to remove the 
doublers on the lattice (see section 7) while B~ always stays positive. Also, only H"*" can become 
gapless and therefore give rise to mismatches in the number of filled states (see section 8). Now we 
shall see that only A'^{U) contributes to the anomaly in the continuum limit: In Figure 10.3 we 
plot A~{U) as a function of j- and see that lim/,^00 •^~{U) = indicating that the state \R— >u 
does not contribute to the anomaly. In Fig (10.4), we plot A'^{U) as a function of j^. Extraction of 
the continuum limit from A{U) yields a value of 0.0254(1) which is in agreement with the righthand 
side of (10.18), namely ^"^ = 0.0254648. Of course, since the anomaly is evaluated as a derivative 
of a functional, it comes in the "consistent" form which for the present abelian computation only 
means that the numerical coefficients in the first row of (10.5) are t- rather than tt-. 

It is rather obvious from the above that if the anomalies cancel between different fermion 
species in the continuum, (as for example they do in the 3-4-5 model - see subsection 10.1) we 
expect the imaginary part on the lattice to be very small for gauge fields typical to small bare 
gauge coupling. We tested this on a 4 X 4 lattice, for the 3-4-5 model, with a set of background link 
variables whose phases were randomly drawn from the interval ( — fo ' 10 ) ^^^ found that the residual 
imaginary part was smaller by a factor of about one hundred than the individual contributions from 
each species. 



46 



10.3 Topological charge 

We explained in section 3 in general terms and in section 8 in detail how the overlap encodes 
the topological properties of continuum gauge field space. The purpose of this subsection is to show 
how this works in simple examples. We shall also show numerically that the lattice definition of 
topological charge is robust under the addition of reasonable amounts of noise to a given smooth 
configuration. This indicates that the classical concept of continuum topological charge when 
latticized by our method maintains its relevance at the quantum level. While this is always assumed, 
it is not a trivial feature because dominating configurations in a path integral are very different 
from smooth continuous fields. 

We start by defining a continuum gauge field on the torus which carries unit topological charge, 
but, unlike an instanton, has uniform field strength. Later, in subsection 10.6, we shall see that 
when the configuration has the topological charge density more localized the needed mismatches 
also occur. Our tests will show that our regularization reproduces the Atiyah-Singer index theorem, 
i.e. the number of effective zero modes depends only on the topology. Even if the field strength 
is uniform the vector potential cannot be a smooth function on the torus because the related 
U{1) bundle is nontrivial. To define the gauge field it is enough to cut the torus along one of 
its main circles, and make the vector potential change discontinuously across the cut by a gauge 
transformation. In the direction that has been left uncut the gauge potential is periodic and hence 
smooth. The gauge transformation at the cut is nontrivial having unit winding when one goes 
around the cut circle. If we now consider the chiral Dirac operator as a mapping from the set 
of fermion fields with a discontinuity across the cut given by the above gauge transformation we 
discover that the images are in the same set. Thus, the mapping itself is entirely smooth in the 
sense that a small deformation in the source field is mapped into a small deformation of the image 
field. When we put this on the lattice the link fields must reflect both the above discontinuity and 
the the fermion boundary conditions.* The following definitions contain a parameter q. Integer 
values of q correspond to topological charge q and smooth Dirac operators, while the intermediate 
values of q correspond to Dirac operators that have real discontinuities on the torus and are thus 
forbidden in the continuum. Their existence on the lattice is the reason that strictly speaking there 
is no topology on the lattice and any configuration of gauge fields can be continuously deformed 
into any other. 

A uniform electric field E = -^ on a torus has a topological charge equal to q. Such an electric 
field is realized on a finite lattice by the following configuration: 

[1 if TCi 7^ X - 1 |. 2'Kniq-. 

Ui{ni,n2)= |^„pr ■ 27rn2q i [{^=1-1' ^^2 (?ii , ?i2 ) = exp [z— ^J (10.26) 

The lattice sites are (711,712) with < ni^2 < L — 1. If q' is an integer then the above configuration 
produces a uniform electric field. The configuration is periodic in the direction 2. The U's living 
on links in the 2 direction implement the linearly raising vector potential needed to generate a 
constant electric flux. When ni goes from X — 1 to zero we note a discontinuity in U2 by a gauge 
transformation. The torus is cut be severing the direction 1 links (X — I,7i2)-(0,7i2). The Ui's 
living on these links implement the boundary gauge transformation on the fermion field. It is easy 
to check that, if q is an integer, the parallel transporters around all elementary plaquettes are equal 



* We remind the reader that we have implicit underlying antiperiodic boundary conditions 
what is being discussed here is on top of those boundary conditions. 
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■2j_q_ 

to each other and to e i^ . Thus, physically, the configuration is smooth. This is not true when q 
is not an integer (consider the plaquette whose corners are (X — 1, X — 1)(-Zv — 1, 0)(0, 0)(0, L — 1)); 
however the configuration is still perfectly acceptable to our construction. The field configuration 
is well defined for any value of q and we would expect the lattice topological charge to change 
discontinuously as one varies q. 
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Figure 10.5 Spectra of H^ as function of q. The integer values of q correspond to smooth gauge 
field configurations of topological charge q. 

In what follows we consider a 6 X 6 lattice and m = 0.9 in H . In Figure 10.5, we plot the 
spectral flow of H~ as a function of q. We see a gap in the spectrum for all values of q and there 
are equal number of eigenvalues, namely 36, on either side of the gap. This is in agreement with 
the proof given in section 8 that H~ has a gap for all gauge fields. Actually, the bounds on the gap 
(|-B| < m) derived there are seen to be very tight.** At q' = the gauge configuration is trivial and 
there are four-fold and eight-fold degeneracies coming from lattice symmetries and spin. These 
degeneracies are lifted when q is turned on. Levels repel when q is varied and the overall range 
stays pretty constant. 

In Figure 10.5, we plot the spectral flow of H"*" as a function of q. Here we see the gap closing 
as one moves away from q = Q. At q' = there are an equal number of positive and negative 



** To be precise one must add finite volume effects that are dependent on the implicit antiperiodic 
boundary conditions. 
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eigenvalues. Between 0.52 < q < 1.55 we have 35 positive eigenvalues and 37 negative eigenvalues. 
In this region both the overlaps u< R — \R-\- >u and u< L — \L-\- >u vanish. The insertion of 
one creation (annihilation) operator in the right (left) handed overlap renders it non-zero showing 
that these gauge configurations indeed have a lattice topological charge of one unit. The plot is 
therefore separated into sectors defined by —2 < q < —1.55, —1.55 < q < —0.52, —0.52 < q < 0.52, 
0.52 < q < 1.55, and 1.55 < q < 2 with lattice topological charges equal to —2, —1, 0, 1 and 2 
respectively. 

We also note that, except for the crossing eigenvalue, the gap is similar in magnitude to the one 
of H~ and both gaps are indeed of the order 2m as expected from the continuum, formal, overlap. 
The range of eigenvalues of H"*" is restricted relative to the one of H~ at the top and bottom. This 
is easily understood as a result of the change of sign in front of m: indeed the difference between 
the H is of order 2m. This squeezes the levels for H"*" into smaller ranges. Note the discernible 
sloping of the gap in Figure 10.5: Extrapolating, it shows that we cannot make q too large and 
maintain the connection to continuum topological charge. This is just a cutoff effect, and we know 
such must exist since there is no nontrivial topology on the lattice. Another example where cutoff 
effects are seen to affect the lattice-topology continuum-topology connection is a configuration that 
has a small concentrated amount of electric flux easily missed by a too coarse embedded lattice. 



X 



o 




Figure 10.6 Same as Figure 10.4, but now at each new value oi q the corresponding link variables 
have been multiplied by independent random phase factors. 
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In a Monte-Carlo simulation the topological charge is a fluctuating quantum variable. There- 
fore topology changing events have to take place during the evolution. The frequency of these events 
determines the relative normalizations of the contributions coming from the different sectors. These 
normalizations are expected to come out in such a way that the equilibrium configurations satisfy 
clustering. Imagine a sequence of configurations labeled by q as in figure 10.5. It is very unlikely 
that such a smooth evolution will ever occur. A slightly more realistic scenario is one in which 
one takes each configuration labeled by a q' (not necessarily an integer) and adds to it noise. We 
took the lattice configurations (10.26), and, for each q multiplied all links by factors e^'^ where 
(j) was randomly drawn from the interval i^ — ^,^). This is a sizable amount of noise but not an 
overwhelming one. The eigenvalue flows obtained from the new set of configurations labeled by q 
are plotted in figure 10.6. The main message is that all the gross features of the noise free flows 
of figure 10.5 survive. In particular, we get all the crossings we had there, thus providing evidence 
for the robustness of our definition of lattice topological charge. 

10.4 Instanton-anti-instanton background 

We continue to look at large but smooth gauge abelian backgrounds in this subsection. We 
now want to see what happens if the total topological charge is zero, but there are lumps of local 
concentration of electrical flux of sufficient size.* As a specific example we consider a single right 
handed fermion in a gauge background made up of an instanton "ring" separated from an anti- 
instanton "ring" along a torus whose other direction is a small circle representing a finite physical 
temperature. We project this configuration onto an Nt X L lattice where 1/Nt is the temperature 
in dimensionless units. The usage of finite temperature is mainly for technical reasons (it reduces 
the size of the matrices we have to diagonalize), but is quite relevant physically also. The link 
variables that define the field configuration are 

Ui{ni,n2) = 1; U2{ni,n2) = \ ^"^"^^^^ if 1 < n2 < n . q < rii < X; < n2 < Nt] 1 < n< L 

{ 1 elsewhere 

(10.27) 

The link variables are uniform in the ^2 direction. All the plaquettes have no electric field except 

for the ones between ni = and ni = 1 and between ni = n and ni = n + 1. Between ni = and 

ni = 1 there is a net flux of 27r and this corresponds to a localized instanton "ring" with charge +1. 

Between ni = n and ni = n+ 1 there is a net flux of — 27r and this corresponds to a localized anti- 

instanton "ring" with charge —1. The total flux is zero and this configuration has no topological 

charge. Continuum arguments will imply that the chiral determinant will decrease as n increases as 

long as n << L: If the distance n (and the size L) were infinite there would be two fermionic zero 

modes and the chiral determinant would vanish. For large but finite distances the absolute value 

of the chiral determinant can be estimated by computing the two smallest eigenvalues of C^C by 

perturbation theory around the two zero modes. These eigenvalues are in turn determined by the 

amount of overlap of the two approximate zero modes localized at the instanton and at the anti- 

instanton. Simply put, there is an attractive potential between an instanton and an anti-instanton 

induced by the chiral fermion. In this particular case of finite temperature, the zero modes decay 

exponentially and the potential is expected to be linear with the constant of proportionality given 



* Throughout the remainder of this section we shall use the terms instanton and anti-instanton 
in a loose sense, meaning a configuration that has a relatively localized lump of topological charge 
density that integrates to unity. Away from the lump the configuration is trivial. 
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by -^. This can be seen by solving the effectively one-dimensional differential equation for the 
zero mode in the continuum. We wish to show that the overlap reproduces the fermion induced 
instanton - anti-instanton interactions correctly. 
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Figure 10.7 Overlap in a background of an instanton "ring" separated spatially by n lattice 
spacings from an anti-instanton "ring" in two dimensions at finite temperature. 

To do this numerically we fix iVy = 4 and m = 0.9 in H^ and compute the determinant 
for the field configuration in (10.27) for a range of lattices with X = 6, 8, 10, 12, 14, 16, 18, 20. On 
each lattice we consider n in the range 1 < n < X/2. The results for the "chiral determinant", 
u< R — \R-\- >u, are plotted in Figure 10.7. For fixed n, we expect the chiral determinant to 
approach a limiting value as X ^ oo. This is seen in Figure 10.7. The limiting values are expected 
to fit an exponential of the form Cexp( — ^) and the curve shown in Figure 10.7 is one such 
exponential obtained by tuning C . The curve smoothly passes through the points for a wide range 
of n showing that the overlap reproduces the expected behavior quantitatively. 
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10.5 Fermion expectation values in an instanton-anti-instanton background 

Although the gauge configuration considered in the previous section has no topological charge, it is 
made up of one localized instanton of charge +1 and another localized instanton of charge —1. As 
mentioned in the previous subsection, formal arguments imply the existence of two approximate 
zero modes one localized at the site of the instanton and another at the site of the anti-instanton. 
This can be uncovered by looking at fermion expectation values in this background. The method 
to compute fermion expectation values was described in section 5. In this case we would like to 
compute the fermion propagator. Normally we would expect the fermion propagators to vanish 
at large distances. But in this gauge background we will get a significant contribution even at 
large distances if the pair of points are located at the sites of the instanton and anti-instanton. In 
particular, from subsection 10.3, we know that there is an excess of negative eigenvalues for H"*" 
in a background with positive topological charge. Therefore we should insert a fermion creation 
operator at the site of the instanton and a fermion annihilation operator at the site of the anti- 



instanton and look at tj< R 



S/9 



\R-\- >u. If we keep y fixed at the site of the instanton field 
(j/ = 0) and vary x, we should see the propagator reach a peak at the site of the anti-instanton. 
Further, this behavior should remain the same, up to an overall multiplication by a factor, if one 
picks another y in the vicinity of 0. 
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Figure 10.8 The fermion two point function in the instanton-anti-instanton background of Fig 
ure 10.7 
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Figure 10.9 Two point function for a = 2,/3 = 1 in Figure 10.8, normalized in cc, plotted for two 
values of y. 



fix L 



We study the above behavior by using the configuration in the previous subsection. Here we 
20 and tc = 10 in addition to fixing Nt = 4. We compute jj< R — \axaa n\R-\- >u- -A-S 
discussed in section 5, we are using Dirac operators to define expectation values of chiral fermions, 
but only half of the modes are physical. This is made very obvious numerically because it turns 
out that for one value of (q:,/3) one gets a matrix element that is significantly larger than the 
matrix elements corresponding to the three other possible (q:,/3) pairs. In Figure 10.8, we plot 



< R 



t 



u ^ It — \axoO-o0\R-\- >(7 as a function of x for all four values of the pairs (a,/3). We see that 
out of the four only a = 2,/3 = 1 contribute significantly. In particular, the complementary pair, 
a = l,/3 = 2 is completely suppressed. Focusing on the dominating matrix element we see that 
indeed the propagator peaks at x = n, the site of the anti-instanton. 

In Figure 10.9, we plot the same propagator in a normalized form, where the square of the 
normalized propagator summed over all x gives unity at fixed y, a,/3. We only plot the a = 2,/3 = 1 
component and display the results for two different values of y, namely y = 1,2, both of which 
are near the site of the instanton. We see that indeed the shapes of the propagators are identical. 
This shows that the fermion propagator in the overlap regularization factorizes into two factors, 
one for the instanton and another for the anti-instanton as expected on the basis of continuum 
manipulations. The factorization is important because it shows how the 't Hooft vertex gets 
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exponentiated in the dilute gas approximation. 

10.6 Fermion expectation values in non-zero topological charge 

The zero modes referred to in the previous two subsections can be computed explicitly in non-zero 
topological charge. Their shape should approximately determine the results we have obtained in 
the previous two subsections and this will provide another check on the selfconsistency of viewing 
the overlap as a discretization of the continuum chiral determinant. In particular, let us focus on a 
gauge configuration consisting of a single localized anti-instanton. The configuration has topological 
charge —1 and can be realized on an Nt X L lattice by the following set of link variables: 

I 1 elsewhere I elsewhere 

The plaquettes between tc^ = and ui = 1 carry identical amounts of electric field adding up to 
a total flux of — 27r. All other plaquettes carry no electric field at all. The Wilson loop encircling 
the torus in the 1 direction at fixed ^2 is different from unity and winds around the complex unit 
circle when ^2 is taken around the torus in the direction 2. The nontrivial winding reflects the 
non-trivial topology of the gauge configuration. Apart from the winding, the configuration is very 
similar to one half of the instanton-anti-instanton configuration considered in the previous two 
subsections. The overlap in such a background is zero due to a mismatch and from subsection 10.3 
we know that H"*" will have one more positive eigenvalue than H~ and therefore the matrix element 
u< R — \axa\R-\- >u should be non-zero. Further, the matrix element, when viewed as a function 
of X can be interpreted as the zero mode in the background gauge field and should therefore peak 
at the site of the anti-instanton. 

In order to establish the existence of the zero mode in the above background we fix Nt = 4, 
L = 20 and expect that one value of a should dominate. Based on the results of the previous 
subsection we know that the relevant a ought to be 2. We plot u< R — \axa\R-\- >u for the two 
values of a as a function of x in Figure 10.10 and see that the a = 2 matrix element is indeed the 
dominant one and that it peaks at the site of the anti-instanton as it should. 

As discussed in the previous subsection, the propagator in the presence of an instanton and 
anti-instanton factorizes and therefore the propagator plotted in Figure 10.8 should represent just 
the zero mode localized at the site of the anti-instanton. Also, the discussions in subsection 10.4 
imply that the decay of the overlap in the presence of the instanton and anti-instanton is governed 
by the decay of the zero mode. The overlap representing the chiral determinant is non-zero only 
due to the nonvanishing overlap between the two zero modes. As one separates the instanton 
and anti-instanton this latter overlap decreases causing the chiral determinant to decrease. It is 
therefore clear that the rate of decay of the chiral determinant should be the same as the rate of 
decay of the zero modes. We can show that these considerations are obeyed by the lattice overlap 
by comparing the normalized zero mode from Figure 10.10 with the normalized propagator from 
Figure 10.9 and with the normalized decay of u< R — \R+ >u from Figure 10.7. The three sets of 
points are plotted together in figure 10.11 and all neatly fall one on top of the other. With this we 
have tested that the real part of the overlap along with the definition of fermion expectation values 
have the right behavior for large but smooth gauge fields in the zero topological sector and also in 
the non-zero topological sectors. 



54 



0.0015 



A 
+ 

8 



V 



0.0010 



t i 



-J — t- 



10 

X 



Figure 10.10 Zero mode in an anti-instanton background. 

10.7 Axial anomaly in nonzero topological sectors 

We have tested the continuum formula for the anomaly in perturbative backgrounds. In this 
subsection we wish to test it in a background that is non-perturbative because it carries topolog- 
ical charge. We cannot do that from the overlap because it vanishes (as does the corresponding 
continuum chiral determinant). However, if we insert a local, gauge invariant operator that carries 
the right amount of anomalous charge we know that its expectation value will become non-zero. 
Suppose we change the fermions by a local gauge transformation corresponding to the anomalous 
(ungauged) group. Since the operator carries anomalous charge the expectation value will change 
even at the classical level in a definite way. At the quantum level, the anomaly will add an addi- 
tional term reflecting the anomaly. We can pick the local inserted operator at a location y say and 
carry out the anomalous gauge transformation in a localized way away from y. Then there is no 
classical variation at all, and all we should get is the anomaly. 

A test of the anomaly in the unit topological sector will support the extension of the Wigner- 
Brillouin phase choice to this case (see section 5.2) and at the same time test the theory in some 
more detail in non-perturbative backgrounds. 

Here we consider a two dimensional abelian vector theory (the matter is a charged Dirac 
fermion) in a background U carrying unit topological charge. Let O be the local operator that has 
a nonzero vacuum matrix element in this background. The overlap now involves the states |F± >u 
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Checking that the zero mode of Figure 10.10 essentially determines the shape of 
Figures 10.7 and 10.9. 



and the definition of their phases is slightly more complicated, as explained in subsection 5.2. The 
expectation value of O in this background is 

i<V-\V->u 



U<V-\V->u 



\U 



<F-|0|F+ >c 



E.„c/<^+|a?„K«] l^+>i 



L lt| 



E.„^<^+KK„]>+>i 



(10.29) 



|F± >u= |-K± >u (^ |iv± >u are the many body states for the vector theory. There are level 
crossings only in the positive side and therefore the modified definition of the phase is needed only 
there. X^^„ a^a \'^xo\ i^ ^^^ '^ Hooft operator that has to be inserted for the overlap to be non- 
zero. It is a Lorentz scalar and is gauge invariant. Since there is an axial anomaly in this theory, 
this expectation value is not expected to be invariant under a local axial phase change made at 
a location far away from O. The anomaly being local is still expected to be given by the gauge 
variation of (10.6), namely, by —i^Ei^x) in the continuum. Further, the middle term in (10.29) 
is invariant under the above axial phase change since O is. (If it were not we would also get the 
classical contribution to the axial variation.) Thus the anomaly must come from the variation of 
the first and third terms. 

As an example of a gauge field configuration with unit topological charge, we take the pertur- 
bative configuration defined in (10.15) with electric field (10.16) and add to that a configuration 
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Figure 10.12 The contribution to the anomaly coming from the "inactive" side, A as a function 
of number of lattice points for a configuration with topological charge. 
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Figure 10.13 The contribution to the anomaly coming from the "inactive" side, A'^ as a function 
of number of lattice points for a configuration with topological charge. 

that has a uniform electric field of strength ^. If we now make an axial gauge variation of the 
type given in (10.19), the uniform field does not contribute to the anomaly. Therefore only the 
perturbative plane wave part contributes and the anomaly should be given by (10.20) with an 
overall factor of 2 to account for two chiralities. To implement this on the lattice, we multiply 
the links representing the perturbative configuration in (10.21) by the ones representing the uni- 
form field configuration in (10.26). We choose the same parameters as in section 10.2 and again 
compute the gauge variations on several lattices. Like in the perturbative case the negative side is 
gauge invariant in the continuum limit and only the positive side has an anomaly. However, the 
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approach to the continuum in the present case is qualitatively different from the perturbative case. 
An extrapolation to the continuum gives 0.0511(2) in good agreement with the continuum value 
^^ = 0.0509296 and the test is successful. 



11. Tests in Four Dimensions. 

In this section we check that the overlap in four dimensions reproduces the abelian anomaly 
and the rules of counting zero modes in abelian and non-abelian backgrounds. We use the transfer 
matrix formalism throughout. For the simple tests we have carried out the computational burden 
is not too bad and since future work with more ambitious goals will most likely focus on the 
Hamiltonian formalism it is useful to establish once and for all that the transfer matrix version 
has been also checked (at least to some extent). The basic logic of the checks is the same as in 
the previous section. We shall see that the overlap formalism passes both perturbative and non- 
perturbative tests successfully. In particular, some of the most crucial ingredients necessary in 
order to reproduce the baryon decay in the minimal standard model are seen to be incorporated in 
the lattice overlap. 

11.1 U{1) anomaly. 

We imagine a single righthanded fermion interacting with an abelian plane wave background. Under 
a gauge transformation the phase of the overlap has to change by the amount dictated by the famous 
triangle diagram. We shall show that the continuum result is reproduced and therefore conclude 
that all the results based on triangle diagrams are reproduced. Combining this with the obviously 
holding Zumino [24] consistency conditions we conclude that all perturbative four dimensional 
anomalies, including non-abelian ones, will be reproduced by the overlap. 

The technical difficulty we have to overcome is that the matrices one has to diagonalize are 
larger in four dimensions. If the system is a lattice torus with equal sides L one deals with matrices 
of size 4X'* X 4X'*. From subsection 10.2 we expect to need to go to X's of the order of 10 at least in 
order to be able to extract the continuum limit. This implies rather large matrices and we do not 
wish to get into ponderous numerics. We opt therefore to use a trick by devising a background that 
has an anomaly but also has translational invariance in two directions, thus reducing the matrices 
down to L^ AL^ X 4X^ blocks. Of course, it is much easier to diagonalize, say, 100 400x400 
matrices rather than one 40000 X 40000 one. We need a configuration that is independent of two 
of the components Ui of the lattice site coordinate n but manages to produce magnetic and electric 
fluxes through a sufficiently large number of plaquettes to make a relatively coarse discretization 
reasonable. The configuration we choose is translational invariant in the 3 and 4 directions: 



TT „«— r COS f ^. TT 



Un,2 = e 



4 

Un,l = e 



COS -Y-{kini + k2n2 + ki)-cos -f (fci rai + fc2'i2) 
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(11.1) 



To get the anomaly we elect to use a more direct approach (also a less efficient one) than the 
one in subsection 10.2. We do this because here, unlike in two dimensions, we shall not carry out 
any other computation that directly tests the imaginary part of the induced action. For this reason 
we put in a free parameter (j) in the configuration. By varying (j), one generates a family of gauge 
equivalent configurations. Using the continuum "consistent" value of the anomaly we conclude that 
we expect the following: 



idT^(U'^,L) 
lim wj 

L^oo TT OCp 



0=0 



-——A3Aikik2 

izTT 



(11.2) 




Figure 11.1 The anomaly as a function of inverse linear lattice size for a four dimensional plane 
wave. 

We set ki = k2 = 1 and A^ = A4 = 0.5 and compute numerically, using the overlap and the 
Wigner-Brillouin phase choice, the left hand side of the above equation for L = 5, 6,..., 12. Anti- 
periodic boundary conditions for the fermions are implicit and the parameter m in the transfer 
matrices has been set to .5. Upon closer examination one realizes that there is more symmetry in 
the problem, related to "rotations" and significantly fewer matrices need to be actually diagonalized. 
Since this is quite technical we omit the details. Also, because of another symmetry, the (j) derivative 
at zero can be computed by evaluating the phase at a single small positive value of (j). 

The results are plotted in Figure 11.1 and the extrapolation to X = 00 is smooth yielding 
— 0.00663(1) for the anomaly in agreement with the continuum value of —0.00663. 

11.2. Topological charge. 

In this subsection we shall discuss a series of tests whose purpose is to ascertain that in four 
dimensions in the presence of nonabelian gauge backgrounds the insertion of the appropriate number 
and kind of creation/annihilation operators between the states making up the overlap renders an 
otherwise vanishing result nonzero. This is done by obtaining the number of eigenstates of the 
"active" transfer matrix whose eigenvalues are less than unity. Actually, the interesting number is 
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generated hy Ts = j [ ^ ^ ) . Since the number of zero modes associated with Ts = 1/2 is 



the deficit or surplus relative to the free case, or, equivalently, relative to the "inactive" transfer 
matrix. Again, the main technical problem is to avoid dealing with too large matrices. 

The gauge group is chosen to be SU{2). Based on our two dimensional experience is is easy 
to guess that on a torus even abelian gauge configurations can carry topological charge. Such a 
configuration can be embedded in one of the f/(l)'s existing inside SU{2), for example the one 
1 ^1 
.0 -1 
equal to the one associated with Ta = —1/2 only even SU{2) topological charge can be generated 

this way. However, in practice we can only look at Ta = 1/2 say, and do not need to carry group 
indices at all. This reduces the size of the matrix by a factor of two. Clearly, it is wasteful to 
fully diagonalize the transfer matrix (we only consider the "active" one, denoted by "+" before); 
it suffices to track eigenvalues crossing unity. We wish to identify all those crossings, but do not 
need the other eigenvalues. In order to be able to focus on the crossings we need a trick to put 
the crossings at the edge of the spectrum of some operator, where they could be observed by 
projecting out an extremal eigenvalue. To achieve this, as explained in section 8, we construct 
T{U) = i[e^ + ('^) + e~^ + ^^> , where the Hamiltonian is the one associated with the transfer 

matrix formulation. Clearly, generic unity crossings in the spectrum of e^+ appear as quadratic 
local minima of the lowest eigenvalue of T{U). The minima of interest occur at one, the lower 
bound of T{U). To obtain the lowest eigenvalue we start with a random unit norm vector and 
operate on it with C — T{U); we renormalize the result and iterate the procedure. C is chosen 
slightly larger than the upper bound of T{U), so that the state we are interested in has the highest 
eigenvalue oi C — T{U). After a few iterations the expectation value of TiJJ) in the state stabilizes 
and this is our estimate for the lowest eigenvalue of T[U) which we denote by j[tQ + j-). This 
method is very unsophisticated but suffices for our needs. 

Using the two dimensional configuration 10.26 in all the 1-2 and all the 3-4 planes we construct 
a four dimensional abelian link configuration that has (in the continuum approximation) fluxes 
Fi2 = -P34 = -^ through every plaquette. q will be varied continuously from zero to 3 and 
whenever it attains an integer value the lattice configuration can be thought of as representing a 
smooth SU{2) configuration with topological charge 2q and uniform topological charge density. 

We set X = 8 and plot in Figure 11.2 the lowest eigenvalue of TiJJ) for several q values 
(represented by plus signs) between zero and three. The eigenvalues are connected by a solid line 
to guide the eye through the flow. One sees clear evidence of three crossings, one in each of the 
intervals connecting two consecutive integers. The number of "zero modes" of the SU{2) fermions 
is twice the one obtained from the figure. The breaks in the line correspond to instances when 
a different eigenstate of e^+, whose eigenvalue t^a minimizes i„ + 7- over n, where n labels all 
eigenvectors of e^+, takes over the dominating position. Since the mapping inverse to i ^ i + j 
is double valued the flow of i + j need not be differentiable, the dominating eigenvector switching 
between the two branches. 

Because q appears in all the exponents of the link variables, and only there, the result for the 
higher q''s can be reinterpreted as evidence that the right number of zero modes will be generated 
for fermions in the representations / = 1/2, 1,3/2 of the gauge group. In particular this includes 
the adjoint and hence the supersymmetric case. It is gratifying to get the right zero mode structure 
in the theory whose continuum limit is supposed to come out supersymmetric because instantons 
play a special role in the dynamics of supersymmetric theories. 

The abelian check is insufficient because we have not tested the simplest case, with unit 
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Figure 11.2 Lowest eigenvalue of T as a function of q for the abelian background. Each time the 
graph touches the value one an eigenvalue of the "active" transfer matrix crosses 
unity. 

topological charge. Also, we have not yet done any test that is sensitive to the nonabelian group 
structure of the transfer matrices. Further, we wish to show that even localized topological density 
will lead to zero modes in the overlap formula. It is inconvenient to work on a torus at this point 
(we shall come back to this later) because the instanton is a solution on four space compactified 
to a sphere. We shall therefore use an open hypercube, i.e. our manifold will look like a piece of 
R^ whose boundary has the topology of S^ and the instanton on the boundary is almost a pure 
gauge. The associated gauge transformation maps the boundary into the group with unit three 
dimensional winding. 

We use antihermitian notation for the continuum gauge fields A^i^x). We wish to latticize the 
following configuration: 



^^^i.^) 



;9\x)d^g[x) g{x) 



X4 



la ■ X 



(11.3) 



In the above equation the ct's are the three Pauli matrices with the indices 1,2,3 combined into 
symbols carrying a vector sign. For R^ sufficiently larger than p^ , the instanton "size" squared, 
the topological charge inside the ball \x\ < R is close to one and it makes sense to expect such 
a configuration to act as if it carried topological charge exactly equal to unity. As a criterion we 
choose to require that R/ p be so large that more than two thirds of the classical instanton action 
in infinite R* be contained in the ball \x\ < R. A simple computations shows that this roughly 
means that we wish the ball radius to be at least three times as large as the instanton size. 

To accommodate the lattice we do not use a ball but rather a hypercube centered around 
X = xq defined by |a; — ccqI < ^ and consider an instanton whose size is 1/3 or smaller. We embed a 
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hypercubic lattice with side L into the continuous medium. We choose L to be even and make the 
center be the point with coordinates {^^, "f'"' "f'"' "f'")' keeping the center off the lattice sites 
at a symmetrical location. The lattice sites x have integer coordinates < cc^ < i — 1. The link 
configuration is computed by evaluating the parallel transporters along the links of the embedded 
lattice associated with the continuum gauge field. We wish to compute the parallel transporters 
exactly. If we do this then a gauge transformed continuum configuration will yield a lattice gauge 
transformed link configuration. Since the eigenvalue problem for the transfer matrices is gauge 
invariant it does not matter which gauge we choose to represent the instanton configuration in and 
we can stick to the "regular" form in (11.3). 

The evaluation of the path ordered exponentials giving the link variables is less formidable 
than one may think because of the symmetry of our configuration. Let us write the gauge field 
configuration more explicitly: 

^.(^) = ;^«.(^)-^ (11-4) 

Denote (cc — xq)^ by A^. Then h{^x) in the above equation is given by h{^x) = /S? -\- p^ and the a's 
are: 



(11.5) 



a^{x) = (A4,A3,-A2), a2{x) = (-A3,A4,Ai), 
a3{x) = (A2, -Ai, A4), a^ix) = (-Ai, -A2, -A3) 

Note that for a given /x the a^(a;) variable does not depend on A^. Since the link variables are 
defined by 

U^{x) = VeIo''^-^''+'^\ (11.6) 

the integrand along the paths does not change its orientation in group space. Only h has some t 
dependence and the path ordering symbol can be dropped. The explicit form of the link variables 
is 

C/^(a;) = e"'''(^)-^*''(^), (11.7) 

where the scalar quantities $^(a;) are given by 



^Jx)= / — = , tan"^ 4 ^^ o (11-8) 



The exponents of the linear combinations of Pauli matrices in (11.7) have familiar closed expressions 
in terms of trigonometric functions. 

In our numerical implementation we wish to work with the smallest volumes possible. It is 
reasonable to demand that we have at least of the order of two lattice spacings sampling the core 
of each instanton and since we want the size of the system to be be at least three times as big we 
take X > 8 (recall that L has to be even and, since the hypercube is open, its linear extent is X — 1 
lattice spacings). We test cases with L = 8, L = 10, L = 12, X = 14 and with p values going from 
— ji for the smallest volume to — g-^- for the largest. 

We do not wish to interpolate the nonabelian configuration to a trivial one this time. Instead 
we vary the mass parameter in the transfer matrix H_|_(C/') as explained in section 8. We need 
to consider only the "active" side, corresponding to negative mass values in our conventions. The 
number of zero modes one would associate with a given gauge configuration and with a given choice 
for the mass parameter m is equal to the number of unity crossings when the mass parameter 
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Figure 11.3 Lowest eigenvalue of the transfer matrix as a function of m for a background consist- 
ing of a nonabelian SU{2) instanton on a hypercube with open boundary conditions. 
Different box sizes L and instanton sizes p are shown. Crossings of the horizontal 
line indicate the creation of an effective zero mode. 



is varied from zero to its final value. This procedure not only tests more directly the method 
described in section 8, but also gives some indication of the dependence of the lattice topological 
charge definition on the extra regulator m. For each intermediate value of the mass parameter 
we construct TiJJ) and project out its lowest eigenvalue in the same manner as above. This time 
however we take the corresponding "eigenvector" and compute the expectation value of e^ + ^^' in 
it. In this way the crossings are made more evident in Figure 11.3. 

We note that one cannot make m too small relative to the instanton size. We also see that tiny 
instantons can easily be missed. The minimal m needed to detect the crossing decreases when we 
increase the fraction of the infinite-space instanton contained in our hypercube. All of these features 
make qualitative sense. Since the open hypercube can be embedded anywhere in a much larger 
system and the associated eigenvectors of the transfer matrix are more or less localized in it, we 
expect that the number of approximate zero modes in a dilute gas of instantons and anti-instantons 
will be correctly reflected by the overlap. 

We still would like to see an instanton-like configuration on the torus performing correctly. 
We also would like to ascertain that the other, higher, eigenvalues of TiJJ) do not come down too 
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Figure 11.4 Several of the lowest eigenvalues of 2T as a function of m for a background consisting 
of a nonabelian SU{2) instanton on a hypercube with toroidal boundary conditions. 
The lattice instanton has been generated by the method of "supercooling". When 
the bottom state hits the value 2 there is a crossing of unity in the spectrum of the 
transfer matrix. 

and mess up the picture of the crossings. We might suspect that they just escaped detection until 
now. To address these concerns we take a configuration generated by the methods of [25] * using 
"supercooling" ("cooling" [26] with the help of a specially augmented lattice action designed to 
stabilize instanton-like configurations). We use anti-periodic boundary conditions for the fermions. 
We carry out almost the same analysis as before, only we substitute the power method used to 
extract the ground state of TiJJ) for each intermediate mass value there by a Lanczos algorithm 
which gives us several low states. We resort back to plotting the estimates for the eigenvalues of 
TiJJ) rather than of H_|_(C/') as we did in Figure 11.3. However, unlike in Figure 11.2 we do not 
include the factor |- in the definition of T[U) so the crossings occur now at 2. The result is shown 
in Figure 11.4 and we see that there is one crossing and that the crossing eigenvalue is clearly 
separated from the rest, as we would have wished. 



* We are grateful to M. G. Perez for supplying us with the configuration. 
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11.3. Shape of zero mode. 

We have stressed before that it is important to see that the shape of the zero modes is repro- 
duced on the lattice. The direct way to do this is to compute matrix elements of the appropriate 
creation or annihilation operators. We did this in two dimensions. A complete diagonalization 
is involved and would require a larger computational effort in four dimensions than we are able 
to invest at present. We opted for a more modest goal: we simply wish to see how the lattice 
states contributing to the lattice topological charge in (8.5) compare with the continuum states 
contributing to the continuum topological charge in (8.6). In the continuum, it is obvious that only 
the zero modes contribute. On the lattice, only the crossing state, evaluated at the crossing point 
contributes. Note that (8.5) holds both for the hamiltonian and for the transfer matrix formulation. 
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Figure 11.5 Shape of crossing state for X = 8, /O = ^ on an open hypercube as a function of 
euclidean lattice distance. The solid lines are ^„2 1 ^■2\?. with C determined by overall 
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In Figure 11.5 we plot the norms squared of the u, v components (the notation is that of section 
8) of the crossing state on an open hypercube of size X = 8 in the background of one instanton of 
size p = J. From Figure 11.3 we see that the crossing takes place at m = —.45. We take the mass 



parameter as m = —0.5. The plot is given as a function of the euclidean lattice distance r = v cc 



and the lines shown are the well known continuum shapes 
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with C determined so that the 
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integral in the continuum over infinite space is equal to the sum over the contributions of both of 
the u, V components over the finite lattice. The ratio of the contributions of the two componets 
u,v on the lattice is preserved between the C-constants corresponding to the two lines. There are 
many lattice points corresponding to a given r^ and all the values are displayed. Their dispersion 
gives a feel for the amount of 0(4) invariance breaking due to the lattice. We see that the zero 
mode is reasonably symmetric and that it falls off with the correct power. We also see that most of 
it is concentrated in the v component, reflecting the definite chirality in the continuum. The shape 
is unchanged as long as the crossing state is well separated energetically from the other eigenstates 
of the "active" transfer matrix. 



12. Discussion 

We claim that our work up to now shows that our regularization reproduces everything we 
know about the continuum field theories we are considering. With this we can declare victory on 
the problem of regularizing chiral gauge theories; all other questions about the theories can be 
answered only after actually carrying out the integrals we have defined. This point of view may 
not be widely accepted and we shall have more to say about this later. First, however, we accept 
the regularization and try to organize our thoughts about the dynamics by listing a few simple 
possibilities. 

Consider an anomaly free gauge theory. The pure gauge action is given by the usual single 
plaquette term. The lattice gauge coupling is /3 and to get to the continuum limit we have to take 
P to positive infinity. The parameter m in the fermionic sector is fixed and does not get tuned. 

The most obvious question is whether there is a phase transition separating strong coupling 
(small /3) from weak coupling. We would like not to have such a transition, at least in the vector-like 
case, but we think that it is more likely that such a transition occurs. The basic reason is related to 
instantons: For large /3 we expect the expectation value of the 't Hooft vertex to be nonvanishing. 
However, for small /3 the link variables fluctuate rapidly and could be approximated by zero; in 
that case, when m is varied from zero to its chosen value, the gap of the "active" Hamiltonian never 
closes and eigenvalue flow across zero does not occur for almost all gauge configurations. Without 
an imbalance in the number of negative energy states the 't Hooft vertex has zero expectation 
value. 

A possible way to investigate this more thoroughly at /3 = is to expand everything in the link 
variables. To leading order the overlap is unity since the two Hamiltonians have identical negative 
energy subspaces. We may therefore work in a "quenched" approximation, in which the effect of the 
overlap on the distribution of the link variables is ignored. (Another way to justify the "quenched 
approximation" is to work in the limit of an infinite number of colors.) We are then interested in 
the properties of the two random Hamiltonians H (JJ). A simple random variable that can exhibit 
transitions is the average level density defined (for the active side say) by p{E, U) = ^- 6{^E — Ei{U)) 
where H^{U)(j)i = Ei{U)(j)i. The simplest observable is the average p{E) =< p(^E,U) >u and a 
possible question is whether p{E) has a gap or does not have a gap around E = Q. One can 
also consider the topological charge observable q{U) = J_ dEe(^E)p(^E, U) where e{E) is the sign 

function. At zero 6 parameter < q{U) >u= but the topological susceptibility, volume" ' ^^^ 
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change discontinuously as /3 varies. The question is whether the topological susceptibility is zero 
or non-zero. 

To summarize, the question is whether a phase transition separating strong from weak coupling 
exists in the quenched approximation and/or in the full dynamical model. The suspected transition 
is somehow related to global anomalous symmetries. In exploratory numerical studies carried out 
in collaboration with Gyan Bhanot [27] we saw some indications from the SU{2) lattice topological 
charge susceptibility that indeed such a transition takes place in the quenched approximation. 

A deeper question relates to strictly chiral theories. If we go back for a moment to an action 
point of view we know that the overlap represents a strictly infinite number of charged fermionic 
lattice degrees of freedom. (In the vector-like case one does not need to be strict about the infinity, 
as long as one tolerates small masses for the fermions.) If we accept that the strict infinity of 
fermionic species is more than an opportunistic trick we may attach importance to the fact that 
Elitzur's famous theorem [28] about absence of spontaneous breakdown of local symmetries does 
not hold necessarily if the number of fields per unit volume is infinite. Do we have an indication 
here that local gauge symmetries can really be spontaneously broken in chiral gauge theories ? 

Of course, before we speak about breaking the gauge symmetry spontaneously we must argue 
how it gets restored on the lattice when the target theory is anomaly free. It is the most important 
dynamical question facing our approach. Before we discuss this point in general terms, we focus 
on a simpler, specific question: For backgrounds that are trivial up to a gauge transformation, 
does one obtain a continuum like behavior in the abelian case, where there is no dependence on 
the gauge in the chiral determinant because the anomaly vanishes for zero field strength. On the 
lattice we might expect some residual gauge breaking terms that vanish only in the continuum 
limit. However, in our formalism the two dimensional gauge invariance on the trivial orbit is exact, 
just like in the continuum. When we try to go to the continuum limit we are sampling mainly the 
neighborhood of the trivial orbit in orbit space, so the above bodes well for the chances to obtain 
the right target theory. We now present the proof of our assertion: 

Let |i2± >ig and |X± >ig be the many body states for a 2-D U(l) gauge field configuration 
which is a gauge transformation of unity and let the states obey the Wigner-Brillouin phase choice. 
Using (7.8) and (7.14) we can write these states as 

\R± >ig -C)^\R± >i ———————— 

\r<R±\Q\R±>r\ 

\L± >ig = g'\L± >i , 1^1 I 

|i< L ± \y\L± >i I 

All the dependence on g comes in through Q. We claim that the overlaps are, in fact, independent 
of g. That is, 

i9<R-\R+>i9=i<R-\R+>i; ig< X- |X+ >ig= i< X- |X+ >i . (12.2) 

To see this, we begin by noting that G is a unitary matrix and from Lemma 4.0, it follows that 

i<ii;±|g|ii;±>i 

The above equation is valid for both ± independently. This is because the full set of single particle 
states that make up |i2± >i and |X± >i form a complete set in which we can represent G. The 
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numerator and denominator are simply the determinant of the two blocks that appear in Lemma 
4.0. The explicit form of the denominator in the above equation is 



i< X ± \g^\L± >i= det^ Y,{ [V'k^(^"; 0) 9*{x)ipf^{xa; 0) 



(12.4) 



There is no group index above because the group under consideration is U(l). We have used the 
notation A = in the right hand side of the above equation to be consistent with the notation in 
section 4. Now we can use (4.24) to rewrite the above equation in terms of the right handed single 
particle states. We get 



i< 



L±\g^\L± >i= det^^{V'f.^(a;a;0)5*(a;)[V'f.±(a;a;0)]* 



(12.5) 



We have used the fact that jf = 1. The phase factors in (4.24) cancel out in the determinant. The 
above equation simply says that 



i< X ± \g^\L± >i= i< i? ± \g^\R± >i 



(12.6) 



(12.6) when used in conjunction with (12.3) says that the phases of i < X ± |^|X± >i and i< 
R ± |^|-ff± >i do not depend on the ± signs and are equal to half of the phase of det G. Therefore 
when we compute the overlaps in (12.2) the phases from the positive and negative sides cancel and 
the result is independent of ^r. 

Having addressed the very special case of the trivial orbit in two dimensional abelian models, 
we now turn to the general picture. The following presents our point of view in more detail, as 
promised in the opening paragraph. 

Our solution to the problem of regularizing anomaly free chiral gauge theories is not perfect 
[29]. A perfect solution would elevate the algebraic conditions for anomaly cancelation to the non- 
perturbative level and show that they are a necessary and sufficient condition for writing down a 
completely gauge invariant well defined path integral. Whether such a perfect non-perturbative 
regularization exists or not is an interesting and potentially deep question.* However, it is logically 
possible that while gauge invariant continuous chiral gauge theories exist, perfect lattice regular- 
izations do not; and even if they do, it is plausible that the desired limits will be approachable from 
imperfect schemes also. 

We claim that our imperfect scheme is quite natural, almost nai've in the sense that it is 
all based on one simple viewpoint, exposed is section 2. Conceptually, it is much simpler than 
other approaches, in particular the ones based on strongly coupled lattice Yukawa models, where 
some non-trivial dynamics seems to be required and the ultimate goal appears to be a perfect 
regularization, although the role of anomalies is not made explicit. 

To argue the simplicity of our approach in more detail we start by reviewing the situation in 
the continuum, at the level of the formal path integral [30]. The integration variables are some set of 
chiral fermions i/j and some gauge fields A. The integrand does not depend on all of the integration 



* In our formalism we are defining the phase of each chiral fermion multiplet separately whereas 
it is sufficient to define the overall phase of the combined many body states. It is conceivable that 
a natural gauge invariant definition can be found when the theory is anomaly free. 
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variables: a gauge transformation ij; ^ i/;^, A ^ A^ leaves the integrand invariant. If we could do 
that in a natural and local way we would eliminate the redundant variables completely and replace 
the set of integration variables by gauge orbits and fermions. We cannot however do that and we 
are stuck with using more fields than we want. The special features of chiral gauge theories can 
be cast into a simple framework if we are willing to accept that the fermion integration measure 
depends on A. Fujikawa has shown that the anomalous dependence of the chiral determinant on 
the coordinates along the orbits in gauge field space can be usefully viewed in this way. Therefore, 
generically, the functional integral is not gauge invariant after all and we do not have a theory that 
describes fermions in interaction with gauge orbits. However, if the anomalies cancel, we do have 
such a theory. 

The clearest way to see the difference is to integrate out the fermions and obtain a functional 
of A, e^^^> . When anomalies cancel T[A) is invariant under A ^ A'^ and when they do not we have 
r(A^) — T[A) = kT^^[A,g) where k is some quantized non-vanishing constant and T^^[A,g) is a 
universal, known, functional of both A and g. While T[A) is non-local and unknown T^^[A,g) is 
both known and local. If anomalies do not cancel we need to deal with one more degree of freedom 
but, essentially, we know the exact dependence of the induced action on it. 

Similarly to the gauge invariant case we cannot explicitly separate the degrees of freedom into 
orbits and coordinates along orbits in a local and natural way. We deal with this problem by adding 
a redundant variable, g, and having the integrand and measure invariant under a (new) gauge 
transformation A -^ A'^, g -^ h~^g. Under this new transformation T{A) + kT^^{A,g){= T{A^)) 
is trivially invariant. We therefore change the bosonic action by adding the term kT^^[A,g) to 
the original action.* We have not really altered the original functional integral because the new 
gauge invariance can be used to set the new integration variable g to unity and then the added term 
vanishes. So up to an uninteresting overall factor, nothing has changed. However, in the new form 
we can deal with the coordinate along gauge orbits first: We can now try to integrate out g because 
we know the exact dependence on it. After that, the new and old gauge invariances mean the 
same thing and we are left with something gauge invariant, i.e. the "true" orbit - orbit interaction. 
Unfortunately, the integration over g changes the theory drastically, and in four dimensions, we 
cannot make sense of the result even perturbatively. Basically, the problem is that g needs to 
be a dimension zero field to honestly play the role of a coordinate along the gauge orbit. In two 
dimensions there is no problem with this and therefore one can construct anomalous gauge theories 
perturbatively, but in four we do not know how to do that, except in the trivial case where there 
is no dependence on g in the path integral at all. In four dimensions, when there are anomalies, 
one can either try to view ^r as a dimension zero field similar to the nonlinear fields appearing 
in effective Lagrangians, but the theory then has a limited energy range of applicability (i.e. the 
theory is non-renormalizable), or, hope that somebody will discover a new way to quantize the 
degrees of freedom g in the T^^[A,g) action [31]. Anyhow, we cannot construct anomalous gauge 
theories around an ordinary gaussian fixed point as we do with all other relativistic field theories 
known to date. 

When we go to the lattice, it is best to avoid having too explicit a Grassmann integration in 
mind, because then we cannot realize Fujikawa's gauge dependent measure [32]. The simplest is 
to start from the chiral determinant itself and require the lattice expression for it to be a faithful 



* The result of integration over g is the same as averaging the overlap over all gauge transforms 
of the link variables U. 
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representation of the continuum e'^^"^) . This is what we have done. When we compute on the lattice 
T[A^) — T[A) we do not get an answer exactly proportional to T^^[A,g) but we know now that 
we do get something sufficiently close to it for slowly varying A and g. For an anomalous theory we 
would need to integrate over g the exponent of the Wess-Zumino action and we do not know what 
the answer would be. Unlike in the continuum, what we have to do is precisely defined, but the 
result, obviously, is still unknown. If the anomalies cancel in the continuum the cancelation on the 
lattice will not be exact. So, we shall still have to carry out the g integral. The action is not special 
now (except being purely imaginary and obeying some discrete symmetries) and the dependence 
on g is weak in the sense that the leading operator has cancelled out. It is therefore reasonable 
to expect the integration over g to simply give the exponent of some generic, local, functional of 
A. By construction, this functional also has to be gauge invariant and if sufficiently local will only 
inconsequentially alter the tree level lattice bosonic action. That lattice gauge (invariant) theories 
are robust under not too large perturbations by gauge breaking terms has been known for a long 
time [33]. Essentially, the lattice, unlike the continuum, has no trouble dealing with dimension zero 
fields: they tend to decouple from the low energy physics if such exists. This is how we expect an 
anomaly free theory to emerge. We have seen that the lattice gauge dependent remnants in two 
dimensions are indeed numerically small for some trial configurations, so it is plausible that the 
breaking is of the inconsequential type. We feel that the real mystery is not that the anomaly free 
theory could come out of this, but what exactly happens to the anomalous theory on the lattice. 
For this we must understand the ^r-dynamics of the lattice Wess-Zumino action constructed via the 
overlap. 

To summarize, we claim that our approach will lead to continuum interacting anomaly free 
chiral gauge theories without any tuning beyond the need to take the bare gauge coupling to zero. 
To be sure, we have not proven this to happen. But we do claim that if indeed chiral anomaly free 
gauge theories exist our claim is very plausible. In addition, one can always add some explicit local 
gauge breaking terms from the beginning and fine tune them. This would be disappointing to us, 
but we should remember that it is done implicitly in ordinary perturbation theory (it cannot be 
done in the anomalous case because the locality requirement cannot be satisfied and without it the 
perturbative scheme breaks down). 

In short, our basic premise is that a Euclidean regularization of a chiral gauge theory is 
successful if it is based on local objects (the operators H^{U) are local in U) and satisfies the 
following criteria: 

1. It is bilinear in the fermions. This means that in a fixed gauge background all fermion corre- 
lations are given in terms of two point functions and zero modes. 

2. It reproduces an acceptable continuum chiral determinant in weak, slowly varying external 
gauge potentials to any number of external legs attached to the fermionic bubble. By "ac- 
ceptable" we mean that universal terms are exactly reproduced and non-universal ones can be 
written as a series in a basis of local operators of increasing irrelevance. Gauge breaking must 
be allowed to occur, but should be restricted strictly to the imaginary part of the induced 
action. Of course, anomalies must be reproduced in their consistent form. 

3. It reproduces instanton effects in some detail, including non-vanishing expectation values for 
't Hooft vertices and correct long distance behavior of fermion propagators. 

4. The regularized chiral determinant is a naturally defined complex function over gauge field 
space that admits generic zeros in the sense of section 4. The regularized chiral determinant 
respects the symmetries under charge conjugation, parity, global gauge transformations, some 
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selected Eulidean 0(4) transformations and chirality interchange as defined in section 4. 

Not all of the requirements must be obeyed, but our construction shows that it is possible to 
satisfy them all. Conspicuously absent from the list is a requirement of exact gauge invariance. If 
our requirements are met we believe that gauge invariance takes care of itself in the anomaly free 
case. 
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